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I.  INTRODUCTION 


The  present  survey  was  undertaken  because  a  need  exists  for  reliable 
and  inexpensive  (light-weight)  airborne  lightning  warning  and  avoidance 
systems.  Two  applications  of  warning  systems  would  be  in  (a)  warnings  of 
distant  storms,  enabling  a  pilot  to  avoid  severe  weather  and  possible 
encounters  with  hail,  ice  and  dangerous  air  motions,  and  in  (b)  warnings 
of  possible  imminent  lightning  strikes  to  the  aircraft  in  electrified 
clouds.  Use  of  the  familiar  airborne  weather  radar  is  relatively  straight¬ 
forward  with  respect  to  (a),  but  not  with  respect  to  (b).  One  disadvantage 
of  airborne  weather  radar  is  its  relatively  high  cost.  Another  is  the 
difficulty  sometimes  encountered,  due  to  attenuation-masking  effects, 
in  identifying  thunderstorms  hidden  behind  heavy  showers. 

Potentially  valuable  adjuncts  to,  and  possible  inexpensive  substitutes 
for,  airborne  weather  radar  are  lightning  detectors.  These  can  be  electri¬ 
cal  or  optical  in  nature.  By  detecting  the  electromagnetic  radiation  emitted 
by  lightning  and  predischarges  ("atmospherics'1  or  more  briefly  "sferics"), 

RF  detectors  can  warn  of  the  existence  of  electrical  activity  at  a  distance. 
The  electrical  activity  is  usually  associated  with  severe  weather.  Distant 
lightning  warning  can  also  be  provided  by  optical  detection  of  lightning 
flashes. 

Another  class  of  detector  provides  near-zone  warnings.  These  are  elec¬ 
trostatic  field  detectors  that  can  be  useful  for  warning  and  avoidance  of 
lightning  strikes  to  the  aircraft.  They  warn  of  the  presence  of  high  elec¬ 
tric  fields.  High  electric  fields  exist,  not  only  in  thunderstorms  that  may 
already  contain  lightning,  but  also  in  non-thundery  clouds  where  their 
presence  may  cause  the  aircraft  to  inadvertently  trigger  a  strike.  It  may 
be  inferred  from  pilot  reports  that  these  triggered  lightning  incidents  out¬ 
number  by  far  the  natural  lightning  strikes. 

Lightning  strike  and  electric  field  statistics,  and  some  physical  back¬ 
ground  on  ligntning  triggering  by  aircraft,  are  considered  in  Sec.  II.  As 
candidates  for  electrostatic  field  detection  instrumentation,  we  consider  in 
Sec.  Ill  field  mills,  radioactive  probes,  and  corona  points.  These  have 
attractive  characteristics,  and  all  have  been  operated  on  aircraft. 
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Section  IV  deals  with  RF  (sferics)  instrumentation  for  detection 
of  distant  lightning.  Among  the  sferics  detectors  considered  as  pos¬ 
sible  candidates  for  adaptation  to  airborne  use  as  bearing  detectors  are 
crossed  loops,  multiple  loops,  a  time-of-arri val  (TOA)  system  and  an 
interferometer  system.  The  crossed  and  multiple  loops  have  already  been 
operated  on  aircraft,  while  the  TOA  and  interferometer  have  only  been 
operated  at  the  ground.  Nevertheless,  the  TOA  and  interferometer  systems 
appear  to  have  high  potential  accuracy  and  could  be  considered  for  adap¬ 
tation  to  airborne  use.  They  need  further  investigation,  however.  The 
multiple  loops  may  also  be  considered,  but  they  also  need  further  inves¬ 
tigation. 

The  crossed  loop  is  a  well-known  direction  finder.  The  commercially 
available  narrowband  version  called  Stormscope  is  apparently  simple  and  low 
in  cost.  However,  field  tests  (both  airborne  and  ground-based)  made  to  date 
raise  questions  regarding  its  accuracy  and  reliability.  Moreover,  details 
of  the  physical  basis  for  its  data-processing  operations  are  not  available. 
Therefore,  more  extensive  tests  should  be  carried  out.  In  this  connection 
the  general  analysis  in  the  present  report  (Sec.  IV  and  the  Appendices)  of 
crossed  loops  as  bearing  and  range  detectors  should  be  helpful  in  recogniz¬ 
ing  errors  and  in  suggesting  improvements. 

Another  commercially  available  crossed  loop  system  (Lightning  Location 
and  Protection)  employs  wide-band  electronics  and  gates  on  the  initial  part 
of  return-stroke  signals.  It  is  ground-based  and  is  presently  deployed  for 
lightning-caused  forest  fire  detection. 

As  a  single-station  range  detector  the  crossed  loop  would  be  expected  to 
be  inaccurate  in  its  present  state  of  development  which  apparently  assumes  a 
1/r  dependence  of  the  magnetic  amplitude.  Inaccuracies  would  be  due  in  part 
to  (a)  lightning  source  variability,  and  in  part  to  (b)  site  error  (skin) 
currents. 

Bearing  errors  analyzed  in  this  report  include  the  effect  of  slant  or 
inclination  in  (non-vertical )  lightning  radiation  channels,  the  effect  of 
ground  reflection  (ionospheric  reflection  effects  are  considered  negligible), 
and  the  effects  of  eddy  currents  induced  in  the  aircraft  skin  that  produce 
distortions  in  the  local  magnetic  intensity  vector.  The  channel  slant  and 
skin  eddy  currents  both  change  the  apparent  source  direction.  Large  errors 
can  be  introduced  when  the  channel  slant  angle  is  large.  (This  type  of  error 
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can  be  minimized  in  the  case  of  return  strokes  by  employing  gated  wide-band 
(GWB)  electronics.)  In  the  case  of  the  other  type  of  error  (site  error)  the 
skin  currents  can  cause  bearing  errors  (site  errors)  of  about  20°  even  if  the 
loop  is  symmetrically  located  on  the  fuselage.  The  error  can  be  larger  if  the 
instrument  is  mounted  near  the  nose,  a  wing  edge  or  a  vertical  fin  edge.  The 
distortion  of  magnetic  field  amplitude  would  cause  errors  in  range  detection 
by  instruments  relating  magnetic  amplitude  to  distance. 

The  narrowband  crossed  loop  also  has  a  180°  ambiguity.  This  can  be 
removed  by  using  an  electric  field  antenna  in  conjunction  with  the  loop.  How¬ 
ever,  electrical  noise  or  uncontrolled  phase  shifts  in  the  electronic  circuits 
may  make  the  system  unreliable.  A  GWB  return-stroke  system  may  not  have  this  need. 

The  multiple  loops,  TOA  and  interferometer  may  be  less  susceptibla  than 
the  crossed  loop  to  bearing  errors  due  to  lightning  slant.  However,  they 
may  be  equally  susceptible  to  site  errors,  and  this  should  be  investigated. 

If  these  instruments  prove  to  have  high  accuracy  and  insusceptibility  to 
errors  with  respect  to  bearing,  multiple-station  versions  of  these  systems 
mounted  on  aircraft  could  be  considered  for  lightning  location  by  triangulation. 

With  respect  to  single-station  range  detection,  4  possible  systems  are 
considered  besides  the  crossed  loops,  namely,  an  electric  amplitude 
detector,  an  H/E  amplitude  ratio  detector,  a  500  kHz  pulse-height-analysis 
detector,  and  a  waveguide-propagation  analyzer  (low-frequency  dispersion 
characteristics  of  the  earth-ionosphere  waveguide).  All  of  these  have 
attractive  features  and  should  be  investigated. 

Optical  systems  as  bearing  detectors  are  considered  in  Sec.  V.  Several 
possible  constructions  are  suggested.  In  principle  an  optical  detector 
would  appear  to  have  a  minimal  bearing  error  for  a  vertical  lightning  chan¬ 
nel,  but  cloud  illuminations  by  reflection/scattering  of  the  light  over 
extended  regions  (e.g.  "heat"  or  "sheet"  lightning)  would  spread  the  appar¬ 
ent  source  and  make  the  direction  correspondingly  uncertain.  An  optical 
system  is  inexpensive  and  simple,  however,  and  can  be  used  for  example  in 
conjunction  with  one  of  the  range  detectors. 

Section  VI  discusses  airborne  radar,  the  masking  problem,  and  relations 
of  precipitation  echoes  to  lightning. 

A  summary  of  all  systems  considered  for  airborne  use  is  given  in  Sec. 

VII ,  together  with  some  tentative  rankings.  The  systems  are  ranked  accord¬ 
ing  to  availability,  simplicity  of  construction,  maintenance,  calibration 
and  interpretation  of  data,  extent  tested,  potential  accuracy,  versatility, 
and  insusceptibility  to  errors.  The  versatility  category  is  of  interest  in 
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that  some  instruments  can  have  more  than  one  function,  e.g.,  both  range 
and  bearing  detection,  or  both  near-range  and  far-range  detection. 

The  appendices  deal  with  "misdirection"  analyses  for  crossed  loops, 
i.e.,  the  bearing  errors  due  to  slant  lightning,  and  due  to  the  distortion 
in  the  local  magnetic  field  by  induced  skin  currents  (causing  site  errors). 

By  a  combination  of  theory  and  experiment  one  may  determine  correction 
factors  for  site  errors  affecting  a  crossed  loop,  for  any  given  airplane 
geometry.  The  results  of  this  determination  would  suggest  optimum  loca¬ 
tions  for  the  placement  of  the  loop  on  the  airplane.  The  correction 
factors  need  be  determined  only  once.  If  this  is  done,  the  site  errors 
can  be  completely  eliminated  by  suitable  adjustment  of  the  electronic 
amplification.  As  an  illustration,  analytical  correction  formulas  are 
given  in  Appendix  B  for  various  simple  analytical  aircraft  models,  includ¬ 
ing  tabulations  for  a  fuselage  modeled  by  a  prolate  spheroid. 

A  better  representation  of  the  geometric  form  of  an  airplane  can  be 
obtained  by  computer  modeling.  An  example  of  this  is  given  at  the  end 
of  Appendix  B,  where  computed  site  errors  for  a  T-39  airplane  are  dis¬ 
cussed.  General  agreement  is  found  between  the  results  of  the  simple  analyt¬ 
ical  model  and  those  from  the  computer  model.  Nevertheless,  more  detailed 
information  associated  with  the  complex  form  of  the  airplane  can  be  obtained 
from  the  computer  model . 
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II.  LIGHTNING  STRIKE  AND  FIELD  STATISTICS 

A.  Strixe  Statistics 

In  the  last  few  years  a  large  amount  of  information  on  lightning 
strikes  to  aircraft  has  been  collected  and  reported  in  the  literature  (Corn, 
1979;  DuBro,  1980;  Fisher  Plumer,  1977).  Newman  and  Robb  (1977)  provide 
the  following  brief  statenv  on  airline  pilot  surveys:  "Recently 
studies  by  Hourihan  (1974)  and  Plumer  and  Perry  (1975)  have  provided 
fairly  extensive  data  on  lightning  strikes  to  modern  jet  airliners. 

According  to  Hourihan,  over  80%  of  the  strikes  occurred  within  clouds, 
in  precipitation  and  in  turbulence.  In  over  half  of  the  strikes,  St. 

Elmo's  fire  (corona)  was  observed  before  the  strike.  Most  strikes 
occurred  near  the  freezing  level,  from  -5°C  to  +5°C,  and  generally  while 
descending  or  ascending  through  about  3  km  altitude."  This  reflects  the 
awareness  of  the  fact  that  lightning  definitely  is  a  hazard  to  commercial 
as  well  as  military  aircraft.  Ways  and  means  have  to  be  found  to  mini¬ 
mize  the  lightning  danger  to  flight  personnel  and  passengers,  the  light¬ 
ning  damage  to  airplanes  and  electronic  equipment.  For  example,  fly-by- 
wire  and  digi tal ly-control 1 ed  aircraft  may  be  extremely  vulnerable.  Prior 
to  evaluating  the  advantages  and  limitations  of  the  different  warning 
systems  to  be  discussed  later,  it  will  be  worthwhile  to  assess  certain 
aspects  of  the  papers  referenced  above  under  the  point  of  view  of  this 
report. 

First  we  may  resolve  a  certain  controversy  between  pilots  who  report 
that  lightning  discharges  are  somehow  caused  by  the  aircraft  and  the 
scientists  who  discount  these  reports  on  the  simplistic  basis  that  suffi¬ 
cient  charge  cannot  be  stored  on  an  aircraft  to  account  for  the  energy 
required.  The  pilot  observations  may  be  found  in  a  United  Airlines  report 
(Harrison,  1965).  The  situation  is  best  described  by  Clifford  (1980)  as 
follows: 
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"The  pilots  almost  unanimously  agreed  that  there  are  two  distinct 
classes  cf  lightning  observed  in  flight.  The  most  common  variety  usually 
occurs  while  flying  in  precipitation  at  temperatures  near  freezing.  This 
type  is  preceded  by  a  buildup  of  static  noise  in  the  communication  gear  and 
the  presence  of  corona  (St.  Elmo's  fire)  can  be  observed  if  the  flight  is 
at  night.  The  buildup  may  continue  for  several  seconds  before  the  discharge 
(lightning  strike?)  occurs.  The  discharge  terminates  the  static  and  corona. 

The  second  variety  occurs  abruptly  with  no  warning.  It  is  most  likely 
to  be  encountered  in  or  near  thunderstorms,  in  contrast  to  the  former  vari¬ 
ety  which  is  more  likely  to  be  experienced  in  precipitation  that  has  no  con¬ 
nection  with  thunderstorms.  Pilots  tend  to  believe  that  the  slow  buildup 
type  of  discharge  is  not  a  true  lightning  strike  but  rather  a  discharge  of 
excess  charge  built  up  on  the  aircraft  by  flying  through  the  precipitation. 
This  non-thunderstorm  type  greatly  outnumbers  the  other.  Both  kinds  can 
create  a  brilliant  flash  and  a  boom  which  can  be  heard  throughout  the  air¬ 
plane. 

The  response  of  scientists  to  the  pilot's  static  discharge  theory  has 
been  universally  negative.  They  insist  that  insufficient  charge  can  be 
stored  on  an  aircraft  to  produce  a  discharge  which  looks  and  sounds  like 
lightning.  Scientists  are  even  more  emphatic  that  insufficient  energy 
could  be  contained  in  such  a  static  charge  buildup  to  produce  any  visible 
evidence  such  as  burn  marks,  pitting  or  other  damage  on  the  aircraft.  Yet, 
the  pilots  continue  to  insist  that  the  aircraft  is  discharging  and  that  the 
discharges  do  manifest  themselves  by  bright  noisy  arcs  and  (not  all  pilots 
are  sure  about  this)  visible  damage.  The  controversy  has  been  character¬ 
ized  as  a  difference  in  view  between  scientists  of  long  standing  and  pilots 
of  long  sitting." 
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We  believe  that  the  aircraft  can  trigger  a  lightning  strike  in  an 
external  field  with  no  cnarge  at  all  on  the  aircraft.  In  support  of  this 
assertion  we  propose  the  following  hypothesis. 

After  the  Apollo  12  incident,  where  the  rocket  was  struck  by  light¬ 
ning  twice,  it  was  pointed  out  by  one  of  us  (HWK)  in  a  letter  to  Kennedy 
Space  Center  personnel  that  the  Apollo  spacecraft  was  not  struck  acciden¬ 
tally  by  two  lightnings  generated  by  a  cloud  that  otherwise  did  not  pro¬ 
duce  lightning  discharges,  but  that  the  rocket  triggered  these  lightnings 
by  itself,  upon  penetrating  high  electric  field  regions  in  the  cloud. 

This  explanation  was  immediately  accepted  by  Kennedy  Space  Center  as  well 
as  the  scientific  community,  because  the  alternative,  namely  that  the 
cloud  waited  just  to  countdown  zero  and  then  fired  two  lightning  bolts  at 
the  rocket  while  not  being  capable  of  producing  other  lightning  discharges 
before  or  after  launch,  is  too  remote  to  be  considered  seriously. 

Thus,  what  the  pilots  call  the  "common  variety"  of  strikes  to  the 
aircraft,  namely  those  that  are  preceded  by  corona  discharge,  are  trig¬ 
gered  lightning  discharges  caused  by  the  airplane  penetrating  regions  of 
high  electric  field.  The  rare  variety,  that  occurs  without  warning,  is 
an  accidental  hit  of  the  airplane  by  a  natural  lightning  that  originates 
somewhere  else  in  a  thunderstorm. 

Since  such  a  distinction  has  far-reaching  consequences  in  dealing 
with  the  aircraft  lightning  problem,  we  will  discuss  briefly  some  basic 
features  of  a  lightning  discharge  that  are  relevant  to  the  triggering 
problem  (but  that  may  not  be  widely  accepted  or  even  recognized).  The 
electric  field  in  a  thunderstorm  is  usually  not  high  enough  to  reach  the 
breakdown  value  of  2000  to  3000  kV/m  (the  lower  values  occurring  at  higher 
altitudes).  However,  a  raindrop,  being  electrostatical ly  similar  to  a 
conductive  sphere,  has  a  field  concentration  factor  of  3  at  the  upper  and 
lower  points  (at  its  "poles")  in  the  direction  of  the  thunderstorm  field. 
Therefore  the  external  field  necessary  to  cause  breakdown  there  is  reduced 
by  a  factor  of  3.  In  addition  to  this  the  raindrop  will  deform  into  a 
spheroidal  shape  under  the  influence  of  the  field,  and  the  field  concen¬ 
tration  factors  at  the  poles  or  tips  will  increase  accordingly.  Values 
of  5  to  10  may  be  reached  by  this  deformation.  This  brings  the  thunderstorm 
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field  required  to  produce  breakdown  at  the  tips  down  to  values  of  300  to 
400  kV/m,  e.g.  values  that  have  been  measured  in  thunderstorms.  If  break¬ 
down  is  reached  at  the  tips  of  an  elongated  raindrop  the  raindrop  will  go 
into  corona  discharge.  In  this  stage  the  raindrop  is  similar  to  an  air¬ 
plane  in  that  corona  discharge  is  observed  by  the  pilots  before  a  light¬ 
ning  is  triggered.  The  corona  discharge  is  a  relatively  stable  discharge 
and  without  a  further  energy  input  the  plane  as  well  as  the  drop  would 
remai-n  in  corona  discharge  until  the  external  field  drops  below  the  break¬ 
down  value.  The  plane  could  obtain  an  additional  energy  input  for  instance 
by  the  impact  of  raindrops  or  ice  crystals  on  the  airplane  skin  or,  in  the 
absence  of  precipitation,  by  entering  a  higher  field  region. 

With  the  additional  energy  input  the  stable  corona  discharge  turns 
into  an  unstable  spark  discharge,  e.g.  a  lightning  discharge.  Note  that 
it  is  not  necessary  that  the  airplane  be  charged  with  an  amount  of  charge 
usually  associated  with  a  lightning  discharge.  A  minute  fraction  thereof 
is  enough  to  turn  corona  discharge  into  a  spark  discharge.  If  this  is 
achieved  the  energy  required  to  let  the  spark  grow  into  a  lightning  dis¬ 
charge  is  supplied  by  the  thunderstorm  field.  This  explanation  should 
remove  the  objections  of  the  scientific  community,  especially  since  the 
same  principle  works  even  for  a  raindrop.  A  raindrop  already  in  corona 
discharge  may  acquire  an  additional  energy  input  by  a  close  approach  to 
another  raindrop,  forming  with  it  a  conductive  filament  of  doubled  lengtn, 
or  by  being  swept  by  turbulence,  updrafts  or  downdrafts  into  regions  of 
higher  field  strength.  If  the  additional  energy  input  is  enough  to  turn 
corona  discharge  into  a  spark  discharge  --  or  in  lightning  terminology 
into  a  streamer  --  this  streamer  has  a  good  chance  to  grow  into  a  light¬ 
ning  discharge.  The  energy  for  this  growth  is  again  furnished  by  the 
thunderstorm  field. 

Needless  to  say,  this  picture  of  the  sequence,  corona  discharge- 
streamer-lightning  discharge,  is  an  oversimplified  sketch  of  a  rather 
complex  problem  in  plasma  physics.  It  is  sufficient  for  our  purpose  here 
to  point  out  the  similarity  of  a  natural  lightning  triggered  by  precipi¬ 
tation  particles  and  the  lightning  triggered  by  an  airplane  labeled  as 
"common  variety"  by  the  pilots.  The  difference  between  these  two  events 
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is  the  greater  capacity  and  field  concentration  factor  of  the  plane.  Con¬ 
sequently  the  airplane  :nay  trigger  lightning  discharges  in  electric  fields 
that  are  muc.’i  ewer  than  those  requ *i red  to  produce  r;a tj r j <  l  > *  .t.  . dis¬ 
charges.  This  fits  in  with  the  pilots'  observation  that  electric  dis¬ 
charges  triggered  by  airplanes  in  non-thundery  clouds  outnumber  those 
occurring  in  thunderstorms.  Probably  hidden  in  this  observation  is  the 
fact  that  thunderstorms  can  be  recognized,  and  therefore  avoided,  by  the 
presence  of  natural  lightning  discharges,  whereas  the  presence  of  high 
electric  fields  cannot  be  detected  by  visual  observation.  If  the  airplane 
goes  into  corona  discharge  (visible  only  at  night)  it  may  already  be  too  late 
for  evasive  action.  In  addition,  there  is  no  information  to  be  gained  from  a 
single  corona  discharge  on  the  fastest  way  to  get  out  of  the  high  field  region. 
It  has  already  been  mentioned  in  an  earlier  report  of  this  contract  that 
all  thunderstorm  warning  devices  based  on  the  electromagnetic  waves  emitted 
by  lightning  discharges  are  blind  to  the  presence  of  electric  fields 
too  low  to  produce  natural  lightning  but  high  enough  to  produce  airplane- 
triggered  lightning  discharges.  Instruments  capable  of  detecting  distant 
lightning  as  well  as  electric  fields  are  the  field  mill,  the  radioactive 
probe  and  the  corona  point,  which  will  be  discussed  in  the  next  section. 

The  variety  of  lightning  strike  that  occurs  abruptly  and  without 
warning  is  rare  according  to  the  pilot  reports.  This  is  the  case  where 
the  plane  is  accidentally  in  the  path  of  a  natural  lightning.  The  light¬ 
ning  is  extremely  shortsighted.  It  would  take  no  notice  of  the  plane  if 
its  intended  path  passes  the  airplane  at  a  distance  greater  than  several 
lengths  of  the  airplane  body,  say  in  rough  figures  greater  than  100  m.  If 

we  define  the  target  area  as  a  circle  with  a  radius  of  r=100  m  =  0.1  km, 

-2  2 

the  target  area  T  is  about  T=3xl0  km  .  Let  us  assume  that  the  cross  sec- 

2  2 

tion  of  a  thunderstorm  where  natural  lightnings  occur  is  S=6x6  km  =  36  km  . 
Therefore,  if  we  assume  that  lightning  is  approximately  equally  distributed 
over  the  thunderstorm  cross-section,  the  probability  P  of  a  hit  is  given  by 
the  ratio  of  target  area  to  thunderstorm  cross-section,  i.e.,  P=T/S=1 / 1 200. 

This  figure  would  be  further  reduced  if  we  take  into  account  that  a  light¬ 
ning  lasts  only  for  about  one  second,  i.e.,  to  produce  a  hit  the  timing  has 
to  be  exactly  right.  This  may  explain  why  a  strike  of  an  airplane  by  a 
natural  lightning  should  be  a  rare  event  indeed. 


-10- 

The  second  observation  of  the  pilots  of  the  abruptness  of'  the  st^ik 
i  os  tne  aata  of  tne  velocity  of  a  leader  strike.  Even  tne  ;;slov.  ste_ 
eacer  aavances  with  a  velocity  of  about  10?  m/s.  Tne  reaction  trr.e  of 
pilot  v  pi  us  the  response  time  of  the  aircraft  to  changes  in  flight  con 
rols;  would  have  to  be  under  1/100  secona  to  recognize  ana  avoic  a 
eader  one  km  away  coming  towards  the  airplane,  assuming  in  addition  tha 
tie  pilot  is  looking  in  the  right  direction  and  can  see  so  far  into  the 
loud. 

In  conclusion  we  may  summarize  the  discussion  aoove  as  follows: 

1.  The  pilot  observations  cannot  only  very  easily  be  recon¬ 
ciled  with  the  scientific  requirements  of  the  mechanism 
of  a  lightning  discharge  including  the  energy  problem, 
but  may  be  considered  as  an  excellent  confirmation  of 
some  hypothetical  assumptions  concerned  with  the  initia¬ 
tion  cf  a  natural  lightning. 

2.  The  emphasis  placed  by  the  pilots  on  the  distinction 
between  a  static  discharge  (triggered  lightning)  and 
the  abrupt  lightning  strike  (natural  lightning)  should 
be  adopted  in  all  statistical  evaluations.  There  are 
vast  differences  between  tne  two  types  of  discharges: 
available  power,  resulting  damage,  number  and  place  of 
occurrence,  and  very  probably  frequency  spectrum.  _rhe 
latter  parameter  may  play  an  important  role  in  the 
assessment  of  vulnerability  of  electronic  equipment. 

3.  The  notion  that  the  airplane  charge  produced  by  impact 
of  precipitation  has  to  furnish  the  charge  for  a  light¬ 
ning  discharge  should  be  abandoned.  However,  the  energy 
furnished  by  this  charging  process  may  play  an  important 
role  in  the  conversion  of  corona  discharge  into  a  streamer 
leading  to  a  triggered  lightning  discharge. 

4.  In  1946,  L.  P.  Harrison  suggested  that  "the  field  distor¬ 
tion  or  augmentation  created  by  the  presence  of  the  air¬ 
craft  may  raise  an  initially  high  but  sub-critical 
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potential  gradient  to  the  level  where  breakdown  occurs 
at  or  near  the  aircraft.  If  conditions  are  suitable 
the  streamer  could  then  continue  to  propagate  between 
charge  centers  and  a  discharge  would  occur.”  (Fitz¬ 
gerald,  1967). 


B.  Field  Statistics 


We  now  turn  to  a  discussion  of  the  field  distribution  in  a  thunder¬ 
cloud  (Figures  1  and  2,  unpublished).  Figure  1  shows  a  graph  of  the 
vertical  gradient  in  a  particular  storm  versus  altitude  or  temperature. 

The  horizontal  axis  gives  the  gradient  from  -200  kV/m  to  +250  kV/m.  The 
vertical  axis  gives  the  altitude  from  about  3.8  to  6  km  and  has  a  temper¬ 
ature  scale  from  +3  to  -12°C.  The  data  were  collected  by  two  airplanes 
equipped  with  field  mills  making  repeated  (26)  storm  penetrations  at  dif¬ 
ferent  altitudes  during  the  NOAA  lightning  suppression  project  (Holitza 
and  Kasemir,  1974;  Kasemir  et  al ,  1976).  The  base  of  the  cloud  was  a  4  km 
altitude.  The  peak  of  the  negative  gradient  of  -150  kV/m  occurred  at 
about  4.9  km  altitude  at  a  temperature  of  -2°C;  the  maximum  positive 
gradient  of  200  kV/n  occurred  at  5.5  km  altitude  or  the  -10°C  temperature 
level.  The  lower  negative  peak  gradient  indicates  the  altitude  level 
where  ground  discharges  would  originate  and  the  upper  positive  gradient 
where  intra-cloud  discharges  would  originate.  However,  both  gradient 
values  were  in  this  case  too  low  to  initiate  lightning  discharges.  It  is 
remarkable  that  during  the  26  storm  penetrations  neither  the  low  altitude 
airplane,  a  T29,  nor  the  high  altitude  B26  triggered  a  lightning  discharge. 
It  is  of  interest  to  note  here  that  the  thicknesses  of  the  high-gradient 
layers  are  rather  small,  of  the  order  of  only  200  to  300  meters.  T\  s 
narrowness  is  even  more  apparent  in  Fig.  2.  The  vertical  axis  again  gives 
the  altitude;  the  horizontal  axis  represents  time  and  is  marked  with  the 
days  from  July  27  to  August  21,  1974.  On  thunderstorm  days  measurements 
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FIG.  1.  Potential  Gradient  vs.  Altitude  in  Thunderstorm 


were  carried  out  of  the  type  discussed  in  connection  with  Fig.  1.  In 
Fig.  2,  however,  only  altitudes  of  cloud  top  and  base  (marked  by  crosses), 
of  maximum  positive  gradient  (marked  by  circles),  and  maximum  negative 
gradient  (marked  by  triangles)  are  given.  The  numbers  above  the  gradient 
symbols  indicate  the  temperature.  The  boundary  lines  of  the  shaded  bands 
mark  the  60  kV/m  values  enclosing  the  maximum  positive  gradient  in  the 
upper  band,  and  the  -60  kV/m  gradient  values  enclosing  the  maximum  nega¬ 
tive  gradient  in  the  lower  band.  We  see  that  maximum  positive  gradients 
have  been  located  between  the  -7  to  -13°C  temperature  levels,  and  maximum 
negative  gradients  between  the  +1  and  -6°C  levels.  The  thickness  of  the 
bands  is  about  1000  feet  =  300  m.  It  is  important  to  note  that  the  tempera¬ 
tures  as  well  as  altitude  levels  of  the  maximum  gradients  are  not  fixed  to 
an  exact  temperature  or  altitude  but  may  change  from  day  to  day.  The 
thunderstorms  were  small  to  moderate  heat  storms  occurring  in  Colorado  and 
Wyoming,  with  ground  level  about  5  kilofeet  above  sea  level.  A  word  of 
caution  should  be  added  regarding  Figs.  1  and  2.  The  blank  area  in  Fig.  2 
between  20  and  28  kilofeet  altitude  does  not  mean  that  this  layer  doesn't 
contain  high  gradients.  It  simply  means  that  this  area  could  not  be  inves¬ 
tigated  because  it  was  above  the  ceiling  of  the  aircraft.  (For  gradient 
data  at  higher  altitudes,  see  Fitzgerald,  1968.)  It  should  also  be  men¬ 
tioned  that  it  is  not  easy  to  locate  the  high  gradient  areas.  Without  the 
records  of  all  three  components  of  the  gradient  provided  by  the  field  mills 
and  by  on-the-spot  evaluation  it  could  not  have  been  done  at  all.  Even  so 
it  is  easily  possible  that  some  high  gradient  areas  could  have  been  missed, 
since  they  seem  to  occur  in  narrowly  restricted  regions.  The  importance  of 
such  data  to  avionics  doesn't  have  to  be  stressed.  The  pursuit  of  this  type 
of  research  to  include  different  storms  (frontal  storms,  severe  storms, 
storms  over  land  with  lower  elevation,  storms  over  sea,  storms  in  different 
seasons,  nonthundery  clouds,  etc.)  is  highly  recommended. 

A  tentative  conclusion  can  be  drawn  from  these  figures  (1  and  2). 

The  thinness  of  the  high  gradient  areas  suggests  that  a  quick  change  in 
altitude  should  be  made  if  an  airplane  finds  itself  in  the  state  of  corona 
discharge.  The  change  should  always  be  to  higher  altitudes  even  though 
there  is  a  chance  of  penetrating  the  higher  positive-gradient  layer.  A 
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triggered  lightning  in  the  lower  negative-gradient  region  may  penetrate 
to  the  ground  and  involve  the  plane  in  a  ground  discharge.  The  lightning 
triggered  in  the  upper  high-gradient  area  will  most  likely  be  an  intra¬ 
cloud  discharge  that  doesn't  have  the  destructive  power  of  a  ground  dis¬ 
charge. 

It  is  interesting  to  compare  these  figures  (1  and  2)  with  the  figures 
"Aircraft  Lightning-Strike  Incidence  vs.  Altitude,"  and  "Lightning  Strikes 
to  Aircraft  as  a  Function  of  Temperature,"  presented  by  Fisher  and  Plumer 
(1977)  and  reproduced  here  as  Figs.  3  and  4.  Most  of  the  lightning  strikes 
in  Fig.  3  occur  in  an  altitude  range  of  5  to  18  kilofeet  and  in  Fig.  4  in 
a  temperature  range  of  -10  to  +6°C.  The  higher  altitude  and  lower  tempera¬ 
ture  boundary  in  Figs.  3  and  4  are  in  excellent  agreement  with  the  boun¬ 
daries  given  in  Figs.  1  and  2.  However,  the  lower  altitude  (5  kilofeet) 
and  higher  temperature  (+6°C)  boundary  in  Figs.  3  and  4  show  a  much  lower 
altitude  than  the  boundary  (14  kilofeet,  +1°C)  given  in  Fig.  2.  This  dif¬ 
ference  can  be  easily  explained  by  the  high  ground  level  of  5600  feet  above 
sea  level  in  Colorado  and  Wyoming.  Therefore  within  reason  the  agreement 
between  the  two  sets  of  data  has  to  be  classified  as  good. 

The  significance  of  the  agreement  between  the  two  sets  of  data  is  in 
the  fact  that  the  occurrence  of  two  physically  different  events  --  lightning 
strikes  to  aircraft  and  high  gradients  within  clouds  --  are  correlated. 

They  are  linked  together  by  the  reasonable  assumption  that  lightning  dis¬ 
charges  are  initiated  or  triggered  in  the  nigh-gradient  areas.  This  result 
suggests  the  possibility  of  conducting  detailed  research  on  the  probability 
of  lightning  strikes  by  using  the  safer  and  more  effective  method  of  locat¬ 
ing  the  high-gradient  areas  in  clouds  rather  than  by  provoking  aircraft 
lightning  incidents  (although  the  latter  is  certainly  a  direct  approach  - 
see  Pitts  et  al  (1979)  and  Maxwell  et  al  (1979)). 

Another  remark  may  be  made  in  this  connection.  Lightning  research, 
testing  of  instruments,  etc.,  are  often  done  in  the  laboratory  with  high- 
voltage  sparks.  This  too  is  a  safe  and  effective  method.  There  is,  how¬ 
ever,  a  marked  difference  between  a  lightning  discharge  and  a  spark.  The 
lightning  is  an  electrodeless  discharge  that  draws  its  power  from  the 
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FIG.  3.  Aircraft  Lightning-strike  incidents  vs  aiutude. 


(Fisher  and  PI umer,  1977) 


to  aircraft  as  a  function  of  temperature. 

and  Plum er,  1977) 
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eiectric  field  of  a  space  charge  distribution  in  the  cloud.  Electro¬ 
technical  ly  speaking,  the  power  supply  represented  by  an  electrified  cloud 
has  a  very  different  internal  impedance  than  a  high  voltage  generator. 

This  internal  impedance  may  exert  a  significant  influence  on  the  electri¬ 
cal  characteristics  of  the  lightning  on  the  one  hand  or  on  the  spark  on 
the  other  hand.  A  similar  difference  exists  in  the  theoretical  modeling 
of  a  lightning  discharge  by  a  transmission  line.  The  transmission  line 
exists  already  before  the  current  pulse  is  injected,  and  resistance, 
capacitance  ar.d  inductance  are  constants  that  can  easily  be  calculated. 

The  lightning  discharge  has  to  build  its  channel  by  itself.  Resistance 
varies  strongly  in  space  and  time  and  is  current-dependent.  Therefore, 
Ohm's  law  (constant  resistance)  doesn't  hold.  In  case  of  the  capacitance 
and  inductance,  even  a  definition  is  difficult.  (Therefore,  use  of  the 
telegraph  equation  in  theoretical  approximations  is  not  sufficient.)  In 
the  theoretical  as  well  as  laboratory  modeling  of  a  lightning  discharge 
it  would  be  highly  desirable  that  some  assurance  be  gained  regarding  how 
far  results  and  conclusions  can  be  applied  to  lightning  problems.  Theory 
would  have  to  be  extended  to  include  plasma  physics,  and  a  check  of  labora¬ 
tory  tests  should  be  made  with  real  lightning  discharges.  This  can  be  done 
by  the  use  of  a  lightning  tower  or  in  our  case  by  the  more  direct  approach 
of  provoking  lightning  strikes  to  test  airplanes  (Pitts  et  al ,  1979;  Max¬ 
well  et  al,  1979). 

Some  relationships  of  lightning  activity  (and  strikes  to  aircraft)  to 
radar  echoes  have  also  been  considered  by  Fitzgerald  (1978).  (See  Sec.  VI 
on  meteorological  measurements.) 


C.  Strike  Statistics  in  Thunderclouds 


Fitzgerald  (1967,  1968,  1970)  has  analyzed  data  from  instrumented  air¬ 
craft  penetrating  thunderstorms  during  the  1965  Rough  Rider  program.  Among 
his  findings  relevant  to  this  section  are  the  following: 
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Very  high  fields  are  not  required  for  strikes  to  aircraft 
(390  kV/m  was  the  maximum  measured). 

The  currents  of  lightning  strikes  to  planes  were  compared 
with  those  to  power  lines: 

Power  lines:  50%  above  15  kA,  10%  above  40  kA 
Planes:  50%  above  2  kA,  10%  above  7  kA 

(largest  current  22  kA,  limited  data) 

In  strikes  to  a  F-100F  aircraft,  in  normal  behavior,  sensors 
indicated  electron  current  moving  upward  from  the  aircraft 
toward  a  positive  charge  region  overhead,  and  electron  currents 
upward  toward  the  aircraft  from  a  negative  charge  area  below 
the  aircraft. 

Most  strikes  occurred  while  the  aircraft  was  in  transition 
from  a  region  with  positive  charge  above  and  negative  charge 
below  into  a  region  with  reversed  polarity  (negative  charge 
above  and  positive  charge  below). 

Strong  fields  (tens  of  kV/m)  are  encountered  at  high  altitudes 
( -68° C  and  50  kilofeet). 

Intense  strikes  were  frequently  associated  with  dissipating 
storms.  Here  we  quote  Fi tzgerald  (1967) : 

"The  data  presented  suggest  that  thunderstorms  in 
their  early  stages  of  dissipation  retain  sufficiently 
large  charge  centers  to  account  for  one  or  more  lightn¬ 
ing  discharges  if  a  suitable  means  of  initiating  a 
streamer  becomes  available.  It  is  likely  that  an  air¬ 
craft  entering  a  storm  in  this  condition  will  act  to 
'trigger'  a  lightning  discharge.  These  clouds  may  have 
little  turbulence  and  no  distinctive  echo  pattern  on  a 
typical  Air  Traffic  Control  radar.  In  normal  IFR  flight 
operations  in  regions  with  thunderstorms  merged  with 
showers  and  cloud  decks,  the  routine  radar  avoidance  of 
the  presently  most  active  storm  portions  may  readily  lead 
to  flight  through  a  decaying  storm  and  the  possibility  of 
an  isolated  lightning  incident  to  the  aircraft." 
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D.  Other  Investigations 

Thunderstorm  penetrations  with  instrumented  aircraft  that  were  struck 
include  those  made  by  Cobb  and  Holitza  (1968),  Nanevicz  et  al  (1977),  and 
Musi!  and  Prodan  (1980).  Their  observations  further  corroborate  those 
considered  in  this  section. 


E.  Corona  and  Triggering 

In  an  interesting  set  of  relevant  experiments,  the  effects  of  corona- 
producing  points  on  the  trigger-breakdown  field  of  the  shuttle-orbiter  were 
investigated  by  Kasemir  and  Perkins  (1978).  They  used  a  scale  model  of  the 
spacecraft  placed  between  the  plates  of  a  large  plate  condenser.  As  part 
of  the  investigation  they  also  used  a  highly-polished  spheroid  to  determine 
the  trigger-breakdown  field  in  the  absence  of,  and  in  the  presence  of, 
corona-producing  points  on  the  spheroid.  Two  principal  results  are  the  fol¬ 
lowing: 

(a)  the  trigger-breakdown  field  is  about  33%  less  with  corona 
points  than  without,  and 

(b)  the  percent  reduction  of  the  trigger-breakdown  field  (33%) 
due  to  the  presence  of  corona  points  is  only  weakly  depend¬ 
ent  on  the  nature  of  the  points  (form,  length,  sharpness,  etc.). 

An  important  ramification  of  these  results  regarding  triggered  light¬ 
ning  strikes  to  aircraft  is  that 

(a)  not  only  is  the  likelihood  of  a  triggered  strike  enhanced 
(33%  lower  trigger  fields)  by  the  presence  of  corona  (in 
addition  to  altitude-dependent  reduction  of  the  trigger 
field) ,  but  that 

(b)  all  aircraft  have  so  many  sharp  metallic  protrusions  acting 
as  possible  corona  points  (antennas,  pitot  tubes,  landing 
gear,  nuts,  lightning  arresters,  exhaust  nozzle  rims,  edges, 
corners,  etc.)  that  they  are  essentially  always  in  corona  in 
strong  fields. 
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III.  FIELD  MEASUREMENTS 

We  consider  in  this  section  3  types  of  instruments  for  measuring 
electrostatic  fields,  all  suitable  for  airborne  use.  We  discuss  first 
the  field  mill,  then  the  radioactive  probe,  and  finally  the  corona-point 
detector,  all  well-known  instruments  that  have  been  used  extensively  for 
many  years.  A  considerable  amount  of  historical  and  technical  data  on 
these  and  many  other  instruments  may  be  found  in  the  comprehensive  texts 
by  Israel  (19/1,  1973)  and  Chalmers  (1967). 


A.  Field  Mills 


The  field  mill  is  an  electrostatic  voltmeter,  of  rugged  design  but 
high  sensitivity.  This  is  an  electrostatic  induction  type  of  instrument. 

A  forerunner  of  the  modern  version  was  designed  by  C.T.R.  Wilson  (1920). 

There  are  many  possible  designs  for  such  an  instrument,  e.g.  planar- 
snutter  field  mill,  cylindrical  field  mill,  rotating  wire,  etc.  In  the 
same  family  of  instruments  as  the  field  mill  are  the  electrostatic  flux- 
meter  (test  plate  moving,  not  fixed),  induction  voltmeter  and  agrimeter 
(test  plate  grounded  when  exposed,  connected  to  measuring  instrument  when 
shielded  -  Chalmers  (1967)).  We  will  discuss  two  types,  the  planar-shutter 
field  mill  and  the  cylindrical  field  mill.  The  instrument  is  capable  of 
measuring  atmospheric  electric  fields  over  a  wide  range,  from  one  V/m  to 
500,000  V/m.  Since  high  electric  fields  are  a  characteristic  feature  of 
thunderstorms,  the  field  mill  is  an  ideal  tool  for  thunderstorm  research 
and  is  widely  used  for  that  purpose.  Its  capability  of  small  and  compact 
design  makes  it  suitable  for  airborne  application. 

Briefly,  in  a  planar-shutter  type  of  field  mill,  a  rotating  grounded 
electrode  alternately  shields  a  fixed  insulated  electrode  from,  and  exposes 
it  to,  the  electric  field  external  to  the  instrument.  A  diagram  of  this 
type  of  field  mill  is  shown  in  Fig.  5a  (from  Chalmers,  1967).  The  a.c. 
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FIG.  5a.  Field  Mill  Schematic  (Planar-Shutter  Type) 
(Chalmers,  1967) 
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voltage  output  of  the  fixed  electrode  is  proportional  to  the  strength  of 
the  field.  The  planar-shutter  type  field  mill  was  designed  to  measure  the 
atmospheric  electric  field  at  the  ground  (see  Israel,  1973).  In  this  case 
the  field  has  only  a  vertical  component.  The  sensitivity  is  about  5  V/m. 

The  instrument  has  a  frequency  range,  the  upper  limit  of  which  is  determined 
by  the  rotation  rate  of  the  shutter  and  the  number  of  vane  sectors,  that  is, 
the  frequency  with  which  the  openings  are  covered  and  uncovered.  Typically, 
for  example,  the  upper  frequency  may  be  10  Hz  for  a  30-Hz  rotation  frequency. 
With  a  capacitive  by-pass,  however,  the  frequency  range  can  be  extended  into 
the  3F  region  and  above  (Smith,  1954).  This  means  that  short  pulses  "field 
changes"  radiated  by  lightning  discharges  can  be  detected.  The  limitation 
here  is  given  by  the  design  of  the  electronic  amplifiers,  and  not  by  the 
sensor  head  (i.e.,  the  field  mill  proper). 

In  general,  the  electric  field  vector  may  be  measured  by  using  multi¬ 
ple  field  mills  with  different  orientations  to  infer  the  three  components 
of  the  field.  At  the  ground,  however,  the  electric  field  vector  is  normal 
to  the  earth's  surface.  Therefore,  one  field  mill  suffices  to  obtain  the 
vector. 

The  extension  of  the  ground  field  mill  instrument  to  airborne  opera¬ 
tion  involves  a  system  of  at  least  4  field  mills,  mounted,  if  possible,  at 
symmetric  locations  on  the  airplane  surface,  e.g.,  the  two  wing-tips,  or 
the  top  and  bottom  of  the  fuselage,  or  similar  positions.  With  an  analog 
or  digital  computer  circuit  the  four  inputs  can  be  analyzed,  yielding  the 
three  components  of  the  external  field  vector,  separated  from  the  field 
generated  by  the  electrostatic  charge  on  the  airplane. 

Airborne  versions  of  the  planar-shutter  field  mill  have  been  designed 
and  flown  by  Gunn  (1933),  Kasemir  (1943),  Clark  (1957),  and  Fitzgerald 
(1965),  and  more  recently  by  C.  R.  Holmes  at  New  Mexico  Tech  and  by  D. 

Olson  at  the  University  of  Minnesota. 

The  cylindrical  field  mill  (Fig.  5b)  was  specifically  designed  for  air¬ 
borne  use  (Kasemir,  1972,  1978).  Its  rotating  electrodes  alternately  shield 
ar.d  expose  one  another,  but  neither  is  grounded.  Its  sensitivity  is  one  V/m, 
which  gives  it  a  larger  range  of  detection  of  thunderstorms  than  the 


FIG.  5b.  Cylindrical  Field  Mill 

Upper  left:  Rotating  split-cylinder  sensor  head  (upper  section), 
on  stationary  base  (lower  section).  White  strips 
are  Teflon  insulators. 

Upper  right:  Opened  sensor  head,  one  cylinder  segment  removed. 

Bottom:  Mounted  on  airplane  nose,  with  eccentric  "humpring" 
charge  compensator  added.  (Dark  line  along  centerline 
is  a  reflection  of  unknown  origin.) 
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planar-shutter  type.  The  basic  frequency  range  is  also  0  to  10  Hz,  with 
the  above-mentioned  extendability  to  RF  frequencies  via  a  capacitive  by¬ 
pass.  The  cylindrical  field  mill  has  an  important  advantage  in  being  much 
less  influenced  by  the  airplane  charge.  This  influence  can  be  reduced 
significantly  by  use  of  a  charge  compensator  or  "hump  ring"  (Kasemir,  1978). 

On  an  aircraft,  two  cylindrical  field  mills  are  sufficient,  one  measur¬ 
ing  two  of  the  three  components  of  the  field  vector,  the  other  measuring 
the  third  component,  plus  one  component  already  obtained  by  the  first  mill. 
This  redundant  measurement  of  one  of  the  field  components  by  both  field 
mills  has  proven  to  be  valuable  and  provides  a  continuous  check  on  the 
proper  operation  of  the  field  mills.  Thus,  two  cylindrical  field  mills  are 
equivalent  to  five  shutter-type  mills. 

The  vertical  cylindrical  field  mill  is  sensitive  to  noise  produced  by 
rain  because  one  of  the  sensitive  cylinders  is  exposed  to  impact  by  rain. 

If  the  shutter- type  field  mills  are  mounted  in  such  a  way  that  none  of  them 
is  facing  the  flight  direction  they  are  protected  from  direct  impact  of 
rain,  and  they  are  comparatively  noise  free. 

An  extensive  planar-shutter  type  field  mill  network  at  the  ground  has 
been  in  operation  at  Kennedy  Space  Center  for  thunderstorm  warnings  during 
the  launch  of  important  rockets,  for  instance  the  Apollo  series,  the  Mariner 
program,  and  so  on.  Many  airplanes  have  been  equipped  with  field  mills  for 
thunderstorm  warning,  and  for  thunderstorm  and  lightning  research  projects 
sponsored  by  federal  agencies.  For  instance,  8  airplanes  equipped  with  field 
mills  participated  at  the  Apollo-Soyuz  launch  at  Kennedy  Space  Center. 

Other  examples  include  projects  conducted  by  the  U.S.  Army  in  1956-1958, 
the  U.S.  Air  Force  Rough  Rider  lightning  research  project,  the  U.S.  Navy 
fog  research  project,  the  NASA  lightning  triggering  project,  the  NOAA 
lightning  suppression  project,  and  the  TRIP  projects  1976-1979.  It  should 
be  mentioned  that  in  the  NASA  lightning  triggering  project  (rockets  fired 
to  trigger  lightning),  in  the  NOAA  lightning  suppression  project  (chaff 
seeding  to  dissipate  the  charge  concentrations),  and  in  the  Air  Force  Rough 
Rider  project  (see  Fitzgerald  (1970),  investigating  probabilities  of  light¬ 
ning  strikes),  the  field-mi  11 -equipped  aircraft  were  used  to  locate  the  high- 
electric-field  areas  in  the  storm  with  the  intention  of  penetrating  and  not 
avoiding  these  areas  of  maximum  lightning  probability. 
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If  an  airplane  approaches  a  thunderstorm  from  a  sufficient  distance, 
the  field  is  essentially  that  of  the  net  charge  of  the  storm.  In  this  case 
it  is  not  difficult  to  display  on  a  cockpit  meter  the  direction  of  the  storm 
using  the  field  vector.  However,  if  the  plane  is  inside  or  close  to  the 
storm  the  field  pattern  becomes  more  complex  and  a  computer  would  be  required 
to  display  the  storm  direction. 

Compared  with  RF  sferics  sensors  (discussed  later.  Sec.  IV),  the 
field  mill  is,  in  its  static  form,  a  short-range  warning  device  capable 
of  sensing  fields  and  thunderstorms  at  distances  up  to  about  10  miles.  At 
this  distance  lightning  discharge  pulses  are  also  clearly  visible  on  the 
record.  However,  using  the  capacitive  bypass  or  RF  mode,  the  d.c.  and  ELF 
part  of  the  signal  is  filtered  out.  This  eliminates  the  capability  of 
recording  thunderstorm  fields  but  increases  the  signal-to-noise  ratio  of 
the  lightning  field.  Therefore,  lightning  discharges  can  be  detected  at 
distances  up  to  100  miles.  It  should  be  mentioned  that  the  use  of  the  d.c. 
and  RF  modes  of  operation  are  not  mutually  exclusive  but  can  be  used  simul¬ 
taneously  in  separate  channels. 

Besides  providing  lightning  information,  the  field-measuring  capabil¬ 
ity  of  a  field  mill  has  two  distinct  advantages.  First,  it  warns  of  the 
possibility  of  triggered  lightning  in  highly-electrified  clouds  that  con¬ 
tain  no  natural  1 ightning (see  discussion  in  Sec.  II).  Such  a  situation 
occurred  in  the  Apollo  12  incident,  during  which  2  lightning  discharges 
were  triggered  by  the  rocket.  Second,  it  warns  of  the  production  of  corona 
discharge  on  the  airplane  that  may  interfere  severely  with  communication. 

The  field-vector  display  in  the  cockpit  would  in  both  cases  indicate  the 
direction  away  from  the  high-field  area. 

An  absolute  calibration  of  an  airplane  field  mill  system  requires  time 
and  effort  at  installation.  Since  the  field  mill  is  used  mostly  as  a  scien¬ 
tific  instrument,  its  calibration,  accuracy,  response,  etc.  problems  are 
fully  recognized  and  dealt  with.  The  same  problems  exist  also  for  the  cali¬ 
bration  of  RF  sensors.  However,  in  this  case  the  necessity  of  calibration 
is  often  disregarded  as  being  of  minor  imoortance,  or  sometimes  not  even 
recognized  as  existing.  This  results  in  bearing  and  range  errors  discussed 
later.  Because  of  the  similarity  of  the  problem  much  can  be  learned  regard¬ 
ing  RF  sensor  calibration  from  field  mill  calibrations.  However,  to  date 
one  has  been  content  here  with  a  coarse  empirical  calibration,  and  errors, 
for  Instance  site  errors  discussed  later,  remain  usually  unrecognized. 


-27- 


8.  Radioactive  Probes 


Radioactive  probes  (collectors)  have  long  been  used  to  measure  the 
potential  of  a  point  in  space  with  respect  to  the  ground  potential,  even 
before  the  electron  tube  was  invented  (see  Israel,  1973).  The  potential 
difference  divided  by  the  height  above  ground  gives  the  potential  gradient, 
provided  the  gradient  is  constant.  The  probe  operates  by  ionizing  the  sur¬ 
rounding  air  within  a  roughly  spherical  region,  of  radius  about  5  cm  at  sea 
level  and  10  cm  at  10  km  altitude.  The  ionized  and  slightly  conductive  air 
discharges  or  charges  the  probe  until  its  potential  is  essentially  equal  to 
that  of  its  near  environment.  Parachute  dropsonde  versions  were  used  in 
1923  and  1940,  and  a  radiosonde  in  1942  (Israel,  1973). 

The  amount  of  radioactive  material  involved  depends  on  the  material  used. 
Typically,  10  microcuries  of  radium  or  100-500  microcuries  of  polonium  may 
be  used.  Polonium  is  a  pure  alpha-emitter  and  is  considered  safe  (used 
commercially  to  prevent  static  charging  on  phonograph  records). 

In  principle  this  is  an  extremely  simple  device.  However,  there  are 
a  number  of  problems  that  have  to  be  considered  in  the  construction  of  a 
field  measuring  device  using  radioactive  probes.  Some  of  these  are  the 
fol lowing. 

The  probe  must  be  well  insulated  from  the  ground  (or  airplane  surface). 
The  insulation  requirement  is  severe:  insulation  resistance  should  not  be 
less  than  10^  to  10  ^  ohm.  The  slightest  film  of  moisture  due  to  humid 
air,  condensation ,  etc.,  would  short-ci rcui t  the  probe  to  ground.  (However, 
Teflon  is  a  hydrophobic  insulator,  and  with  it  the  problem  may  not  be  too 
severe  (R.  Markson,  personal  communication).) 

The  load  placed  on  the  radioactive  probe  must  have  impedance  of  the 
12 

order  of  at  least  10  ohm  so  that  the  probe  can  maintain  its  open-circuit 
voltage.  In  presently-avai lable  technology  this  dictates  the  use  of  field- 
effect  transistors  (FET's)  as  first-stage  amplifiers.  These  amplifiers  can 
handle  comfortably  input  voltages  in  the  range  from  millivolts  to  several 
tens  of  volts.  However,  they  would  be  damaged  by  50  V  or  100  V  at  their 
input  terminals.  On  the  other  hand,  a  radioactive  probe  located  at  50  cm 
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above  ground  would  supply  in  normal  fair  weather  of  100  V/m  an  input  volt¬ 
age  of  50  V.  In  a  moderate  thunderstorm  field  of  10  kV/m,  it  would  supply 
an  input  voltage  of  5  kV. 

Hence  the  FET  amplifier  requires  a  high-ohmic  voltage  divider  to 
reduce  the  input  voltage.  However,  high-ohmic  resistors  have  not  only  a 
temperature  coefficient  but  also  a  voltage  coefficient.  That  is,  their 
resistance  value  drops  significantly  for  voltages  higher  than  about  1000  V. 

In  addition,  the  radioactive  probe  itself  saturates  at  fields  on  the 
order  of  5  kV/m. 

In  spite  of  all  these  problems  the  radioactive  probe  is  used  as  a 
field  meter  since  it  also  has  some  unique  advantages.  It  and  the  corona 
detector  (below)  are  to  our  knowledge  the  only  field-measuring  devices 
that  have  no  moving  parts.  It  is  of  light  weight,  and  has  negligible  power 
requirements  including  the  electronics.  Its  construction  costs  are  very 
moderate.  It  is  essentially  a  fair-weather  instrument,  but  its  capability 
for  measuring  small  fields  down  to  one  V/m  or  less  would  make  it  useful 
for  detection  of  thunderstorm  fields  at  a  large  distance. 

Radioactive  probes  have  been  operated  successfully  on  aircraft  out¬ 
side  clouds  by  Vonnegut  et  al  (1961)  and  Markson  (1976).  In  the  only  air¬ 
borne  comparison  test  known  to  us  in  which  a  radioactive  probe  system  and 
a  field  mill  system  were  flown  simultaneously  on  two  separate  airplanes 
flying  side  by  side,  the  radioactive  probe  data  and  field  mill  data  tracked 
each  other  extremely  well  (Kasemir,  1978). 


C.  Corona  Detector 


The  relation  between  the  current  I  flowing  from  a  corona  point  and  the 
electric  field  E  producing  the  corona  discharge  is  given,  according  to 
Whipple  and  Scrase  (1936),  by 

I  -  a  (E2  -  Eq2);  E  >  Eq  (1) 


-  t  *. 


where  a  is  called  the  "corona  constant"  and  has  an  empirical  value  in  the 
11  -16  2  2 

range  10  to  10  A-m  /V  .  Corona  discharge  does  not  start  until  the 
field  at  the  corona  point  exceeds  the  breakdown  value  in  air.  The  ratio 
of  the  field  at  the  point  to  the  corona-producing  field  E  is  the  field 
concentration  factor.  This  factor  is  fairly  constant  for  the  individual 
point  but  it  depends  on  the  height  of  the  point  above  ground,  usually 
given  by  the  length  of  the  rod  carrying  the  point,  and  on  the  radius  of 
curvature  at  the  point.  Therefore,  large  differences  in  the  field  con¬ 
centration  factor  exist  from  point  to  point. 

The  field  concentration  factor  determines  the  onset  field  E  ,  which 

o 

has  a  value  of  about  one  kV/m  for  a  field  concentration  factor  of  3000. 

If  the  external  field  E  is  smaller  than  the  onset  field  the  corona  cur¬ 
rent  is  zero  (not  negative).  The  onset  field  and  corona  constant  for 
positive  fields  differ  by  about  10  percent  from  those  for  negative  fields. 
Moreover,  these  parameters  depend  on  the  wind.  This  wind  dependence  in 
airborne  applications  has  been  thoroughly  investigated  by  Chapman  (1958). 

An  intensive  study  of  corona  discharges  has  been  made  by  Loeb  and 
his  co-workers  (Loeb,  1965). 

Balloon-borne  corona  points  were  used  as  early  as  1937  and  1941  by 
Simpson  and  Scrase  (1937)  and  Simpson  and  Robinson  (1941)  to  determine  the 
polarity  of  electric  fields  in  thunderclouds.  More  recently  measurements 
of  thunderstorm  fields  using  corona  points  were  made  by  Weber  and  Few 
(1978)  using  balloon-borne  ("coronasonde")  payloads,  and  by  L.  Ruhnke 
(1971b)  using  Mighty  Mouse  rocket-borne  payloads.  At  Ruhnke' s  suggestion, 
R.  Markson  (personal  communication)  complemented  his  airborne  radioactive 
probe  system  by  an  equivalent  corona-point  system,  in  order  to  extend  the 
field-measuring  range  of  the  radioactive  probes  to  higher  field  values. 
Markson 's  results  showed  that  the  two  systems  tracked  one  another  well  in 
their  overlap  range  (R.  Markson,  personal  communication). 

The  combination  of  these  two  devices  (radioactive  probes  and  corona 
points)  deserves  further  investigation.  Both  systems  are  simple  and 
light-weight.  The  corona  points  have  the  advantage  of  being  operated  in 
the  short-circuit  mode  and  therefore  do  not  need  high-ohmic  insulation. 
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They  have  the  additional  advantage  of  being  rugged  c  d  fairly  immune  to 
rain  or  humidity.  Furthermore  they  would  be  more  sensitive  to  corona 
discharge  than  the  communication  equipment  of  the  airplane,  and  can  be 
designed  as  a  system  to  measure  all  three  components  of  the  electric 
field  vector. 


D.  Extension  to  Long-Distance  (Sferics)  Detection 

As  with  field  mills,  radioactive  probes  and  corona  detectors  are 
primarily  short-range  warning  devices,  but  can  be  extended  to  also  detect 
RF  pulses  from  distant  lightning.  The  extension  again  consists  of  using 
capacitive  by-passes.  The  rod  used  to  mount  the  radioactive  probe  or  the 
corona  point  would  serve  as  an  antenna  for  lightning  RF  signals. 
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IV.  SFERICS  DETECTORS 

The  electromagnetic  radiation  emitted  by  a  lightning  discharge  is 
usually  known  to  radio  or  TV  users  as  "atmospherics"  or  "sferics".  It 
can  be  used  in  different  ways  to  determine  the  direction  of  the  source  of 
the  emission,  i.e.,  the  lightning  discharge.  We  will  consider  here  crossed 
loops.  TOA  systems  and  interferometer  systems  will  be  discussed  in  a  later 
section  ("Special  Systems"). 

A.  Crossed  Loops  as  Direction  Finders 

Low-frequency  direction-finding  by  narrowband  magnetic  loops  is  a 
rather  old  and  well-known  technique  (e.g.  Watson  Watt  and  Herd,  1926).  The 
theory  of  the  magnetic  crossed  loop  is  as  follows. 

A  time-varying  magnetic  field  penetrating  the  area  enclosed  by  a  wire 
loop  will  induce  in  the  wire  a  voltage  and  current  proportional  to  the  time 
variation  and  strength  of  the  magnetic  field  and  to  the  area  of  the  loop. 

If  the  normal  vector  of  the  loop  surface  and  the  magnetic  field  vector  are 
parallel  the  induced  voltage  is  a  maximum.  If  the  two  vectors  are  orthogonal 
the  induced  voltage  is  zero.  At  an  arbitrary  angle  of  incidence  a  the 
induced  voltage  depends  on  the  cosine  of  the  angle  a.  If  there  are  two  loops 
whose  planes  are  at  right  angles  to  one  another  (crossed  loops)  the  induced 
voltages  in  the  two  loops  are  proportional  to  cosa  for  one  loop  and  to  siria 
for  the  other  loop.  The  ratio  of  the  two  loop  voltages  is  a  function  of 
the  incidence  angle  a,  but  not  of  the  magnetic  field  strength  since  this 
cancels  out.  From  this  ratio,  therefore,  the  azimuthal  direction  of  the 
lightning  discharge  is  easily  determined.  Note  that  only  two  components 
of  the  magnetic  vector  can  be  detected  with  a  crossed-loop  antenna.  Since 
the  antenna  is  usually  mounted  on  top  of  the  fuselage  the  visibility  of  the 
crossed  loop  is  limited  to  the  horizontal  plane.  Only  D.  Kohl  (1966)  has 
attempted  to  obtain  all  3  components  of  the  magnetic  vector  with  his  n.jlvi- 
loop  antenna.  (See  discussion  on  p.  41.)  There  remains,  however,  a  180- 
degree  ambiguity  in  the  two-loop  system.  This  can  be  resolved  by  using 
an  omnidirectional  electric  field  antenna  in  conjunction  with  the  crossed 
loops.  However,  electrical  noise  or  uncontrolled  phase  shifts  in  the  elec¬ 
tronic  circuits  may  make  the  system  unreliable. 
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1.  Use  of  Electric  Antenna  to  Resolve  180-degree  Ambiguity 


A  simple  way  of  seeing  how  correlation  of  the  phases  of  the  electric 
and  magnetic  vectors  can  be  used  to  resolve  the  ambiguity  is  the  following. 


Consider  one  loop  of  a  narrowband  crossed-loop  system  aligned  with  two 
sources,  one  in  "front"  and  one  in  "back".  In  Fig.  6,  the  plane  of  the  loop 
is  perpendicular  to  the  plane  of  the  paper  and  is  aligned  with  the  sources 
(upper  dot  =  front,  lower  dot  =  back).  The  sources  are  assumed  to  be 
electric  dipoles  in  the  same  planp  as  the  detector,  but  oriented  perpen¬ 
dicular  to  the  plane  of  the  paper'.  Four  sets  of  vectors  are  shown,  at 
(a),  (b),  (c)  and  (d),  with  E,  H,  and  P  in  each  set  denoting,  respectively, 
the  electric  vector,  magnetic  vector,  and  di rection-of-propagation  vector, 
respectively.  Sets  (a)  and  (b)  correspond  to  radiation  from  the  front, 
while  sets  (c)  and  (d)  correspond  to  radiation  from  the  back.  Let  signs 
be  associated  with  E  and  H,  such  that  E  pointing  upward  and  downward  denotes 
a  positive  and  negative  amplitude,  respectively.  While  H  pointing  right  and 
left  denotes  also  a  positive  and  negative  amplitude,  respectively  (readily 
convertible  into  loop  voltage  signs).  The  four  possibilities  may  be  tabu¬ 
lated  as  follows. 


a) 

b) 

c) 

d) 


E_  H_  Source 

+  front 

+  -  front 

back 

+  +  back 


It  is  evident  from  the  table  that  E  and  H  have  opposite  signs  if  the  source 
is  in  front,  while  they  have  the  same  signs  if  the  source  is  in  back.  It 
js  not  difficult  to  apply  the  same  procedure  to  the  second  loop. 

Thus  it  is  readily  seen  why  electrical  noise  can  produce  180°  errors. 
If  the  signal-to-noise  ratio  is  too  low  or  the  electronics  is  ineffective, 
electrical  noise  pulses  of  the  wrong  sign  can  overwhelm  the  electric  field 
amplitude  and  make  it  appear,  say,  negative  when  it  should  be  positive,  and 
vice  versa.  Similar  errors  may  be  caused  by  unintentional  phase  shifts  in 
the  signal  amplifiers. 


/ 
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(a) 

P 


P. 

(c)  A 

H  ®E 


•  source  in  front 


•  source  in  back 


FIG.  6.  Resolution  of  18Q-degree  ambiguity 


2.  Errors  in  Direction-Finding 


In  theoretical  discussions  it  is  usually  assumed  that  the  electro¬ 
magnetic  field  is  linearly  polarized,  implying  that  the  source  current 
flows  in  a  straight  line.  The  lightning  channel,  however,  is  not  a 
straight  line  but  generally  follows  a  rather  tortuous  path  with  branches 
in  all  directions.  Such  a  source  current  has  many  significant  horizontal 
components.  Consequently  the  assumption  of  linear  polarization  can  result 
in  rather  large  errors  in  direction  determination. 


(a)  Slanted  Lightning  Channel  Errors 

Even  if  the  current  can  be  assumed  to  flow  in  a  straight  line,  in 
general  it  is  nonvertical ly  oriented.  If  we  are  restricted  to  detection 
of  two  magnetic  field  components  (parallel  to  the  ground  plane),  the 
"slant"  lightning  channel,  if  it  is  at  any  significant  altitude,  plus  its 
image  in  the  conducting  ground  plane,  will  cause  the  crossed  loop  to  deter¬ 
mine  an  erroneous  direction.  That  is,  the  magnetic  vector  at  the  crossed 
loop,  instead  of  being  perpendicular  to  the  line  of  sight  (as  in  Fig.  6) 
will  have  another  direction  and  a  different  source  direction  (perpendicular 
to  the  observed  magnetic  vector)  will  be  inferred. 

This  effect  has  been  analyzed  by  Kalakowsky  and  Lewis  (1966)  and  Uman 
et  al  (1980),  but  only  with  the  crossed  loop  detector  at  the  ground.  This 
results  in  simplifications.  The  more  general  case  where  the  detector  is 
also  at  a  non-negl igible  altitude  (i.e.  airborne)  above  a  conducting  ground 
plane  does  not  appear  to  have  been  considered  in  the  literature.  Hence  we 
present  an  appropriate  analysis  in  Appendix  A.  In  the  following  we  discuss 
some  key  results  and  sample  applications  of  the  analysis. 

The  error  in  direction  due  to  a  slanted  lightning  channel  may  be 
expressed  as  (see  Fig.  A-l ,  App.  A): 


G  P/P 

Misdirection  Angle  =  arc  tan  y yp1 


(2) 
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where  P.,  P  ,  and  P  are  the  components  ot  the  lightning  dipole  radiator 

x  y  z 

(P  horizontal  in  direction  of  observer,  P  horizontal  and  transverse  to 

X  y 

this,  Pz  vertical),  and  where  G  is  a  geometrical  factor 

G  =  (H  +  Qh)/D  (3) 

where  0  is  the  horizontal  distance  from  the  dipole  to  the  detector,  H  is 
the  altitude  of  the  dipole,  and  h  is  the  altitude  of  the  aircraft,  both 
above  a  conducting  ground  plane.  Q  is  given  by  Eq.  (16)  of  App.  A,  re¬ 
expressed  here  as  follows: 

Q  ■  [(xj  -  x^)P  +  (x?  -  x|)P/c]/C(x*  +  x2)P  +  (x?  +  x|)P/c]  (4) 

where  x-j  =  [D2  +  (H-h)2]  ,  x2  =  [D  +  (H+h)  ]  and  P  and  P  denote 

first  and  second  time-derivatives  of  the  moment  P.  Note  that  Q  vanishes 
when  the  detector  is  on  the  ground  plane  (h=0)  and  G=H/D  from  (3)  then 
leads  to  the  simple  result  given  by  the  cited  papers. 

As  an  example  we  assume  that  the  lightning  is  at  D=15  km  from  the 
aircraft  and  at  height  H=3  km.  We  also  assume  that  the  dipole  is  tilted  at 
45  degrees  such  that  P^=PZ  and  Px=0.  If  the  aircraft  is  on  the  ground 
we  have  h=0  so  that  G=0.2.  Then  from  (2)  the  misdirection  angle  is  arc  tan 
G  =  arc  tan  0.2  =  11.3  degrees.  Next  assume  the  aircraft  is  at  3  km  alti¬ 
tude  (  =  h).  We  may  evaluate  Q  from  (4)  if  P(t)  is  known. 

Or,  we  may  assume  that  the  dipole  is  a  harmonic  radiator  of  fixed  fre¬ 
quency  and  use  the  alternative  complex  form  (24)  of  App.  A.  Thus  we  obtain: 


Re  Q  = 

0  for  h  =  0 

(5) 

Re  Q  - 

(3  +  2k2  D2)  f0r  h  =  H  «  D 

(1  +  k2  D2)  D2 

(6) 

Re  Q  = 

.  (°.81  *  0-8|  02)  for  „  .  „  ,  0 
(1  +  0.59  k2  D2) 

(7) 

(Re  =  real  part) 

where  k  is  the  wave  number.  Note  from  (5)  -  (7)  that  the  correction 
factor  Re Q  for  finite  aircraft  height  (h)  has  negative  values  ranging 
from  zero  to  about  -3,  for  all  values  of  kD.  Here  we  are  assuming  0 


greater  than  H  (distance  from  lightning  greater  than  lightning  height). 
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Now  assume  a  frequency  of  50  kHz  (Stormscope  central  frequency) .  With 
the  lightning  at  distance  15  km  (=0),  at  height  3  km  (=H),  and  with  the 
aircraft  at  altitude  3  km  (=h),  we  use  (6)  to  obtain  Re  Q  =  -  .10,  a  rela¬ 
tively  small  correction  leading  to  G=  0 . 1 8  as  opposed  to  the  h=0  value  of 
0.2.  Hence  the  misdirection  angle  is  nearly  the  same  (under  the  assumed 
conditions)  when  the  aircraft  is  at  3  km  as  when  it  is  on  the  ground. 


(b)  Site  Error 

Another  kind  of  error  leading  to  errors  in  direction-finding  by  mag¬ 
netic  crossed  loops  is  re-radiation  (or  scattering)  of  the  incident  magnetic 
field  by  nearby  conducting  surfaces  (site  error).  When  the  crossed  loop  is 
mounted  on  an  aircraft  conducting  surface  this  error  can  be  important.  We 
are  unaware  of  any  treatment  of  this  effect  for  airborne  application 
(although  the  work  of  Horner  (e.g.  1954)  is  well  known  for  treating  site 
errors  at  the  ground  due  to  buried  cables,  hills,  etc.).  An  analysis  of 
this  effect  for  airborne  application  is  given  in  Appendix  B.  We  believe 
that  this  analysis  is  original.  Some  key  results  and  sample  applications 
follow. 

The  misdirection  error  due  to  the  reradiation  from  the  aircraft  skin 
depends  on  the  geometry  of  the  aircraft  and  the  position  where  the  detector 
is  mounted. 

In  Appendix  B  we  have  assumed  the  magnetostatic  limit  (wavelength  >> 
aircraft  dimensions)  and  solved  Laplace's  equation  to  obtain  the  perturba¬ 
tion  field.  For  an  arbitrary  angle  of  incident  direction  the  longitudinal 
and  transverse  magnetic  field  components  are  enhanced  by  different  factors 
giving  rise  to  misdirection  or  site  error. 

Assume  the  detector  is  mounted  on  top  of  a  fuselage,  and  that  the 
fuselage  has  the  shape  of  a  long  cylinder  of  circular  cross-section. 

Assume  the  incident  direction  is  at  0°  if  along  the  axis,  and  at  90° 
if  perpendicular  to  the  axis.  Then  Table  1  (App.  B)  shows  that  the  error 
is  zero  degrees  if  the  incident  direction  is  0°  or  90°,  but  increases  to  a 
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maximum  of  -19.5°  angle  error  at  an  incident  direction  of  55°.  This  result 
is  independent  of  the  cylinder  radius. 

Now  assume  that  the  fuselage  is  modeled  by  a  prolate  spheroid  of  aspect 
ratio  0.1  (ratio  of  minor  to  major  axes).  This  approximates  an  F-106  Delta 
Dart,  for  example.  Table  3a  (App.  8)  shows  that  a  detector  on  the  top  at  the 
center  of  the  fuselage  experiences  a  maximum  angle  error  of  -18.4°  at  inci¬ 
dent  direction  55°.  If  the  detector  is  moved  forward  (or  backward)  to  very 
near  the  nose  or  tail  (Case  y=0.99b  in  Table  3a)  the  angle  error  increases 
to  -23.6°  at  incident  direction  60°,  and  reaches  -60°  (that  is  the  negative 
of  the  incident  direction)  on  the  nose  or  tail. 

If  the  aspect  ratio  is  0.2  (more  nearly  characteristic  of  a  C130  geom¬ 
etry  which  lies  between  0.1  and  0.2),  then  Table  3b  shows  that  the  maximum 
error  angle  is  -16.4°  at  incident  direction  55°,  if  the  detector  is  centered 
front  to  back.  This  error  increases  to  -30.7°  if  the  detector  is  nearly  on 
the  nose  or  tail,  at  incident  direction  60°.  Again  the  error  is  minus  the 
incident  direction  on  the  nose  or  tail. 

As  the  fuselage  becomes  thicker  compared  with  its  length  the  error 
becomes  smaller  at  the  central  position,  and  larger  at  the  far  forward 
and  aft  positions.  The  model  becomes  a  sphere  when  the  aspect  ratio  is 
unity  (Table  3d) . 

Appendix  B  presents  formulas  for  the  above  errors.  Also  treated  is 
the  case  of  a  cylinder  of  elliptic  cross-section ,  which  can  be  used  to 
model  the  situation  where  the  detector  is  on  a  symmetric  wing  or  on  a 
tail  fin.  When  the  detector  is  on  the  center  of  the  flat  part  of  the  wing, 
away  from  edges,  the  error  is  negligible.  However,  if  it  were  mounted  on 
top  of  a  thin  body  such  as  the  edge  of  a  tail  fin  the  maximum  error  becomes 
very  large. 


(c)  Ionospheric  Reflection  Error  at  Very  Low  Frequencies 

When  the  lightning  source  is  sufficiently  far  away  the  radiation  can 
propagate  by  multiple  paths.  One  is  a  ground  wave  which  arrives  first. 


This  is  followed  by  the  first  ionospheric-reflected  wave  ("sky-wave")  and 
then  by  many  more  echoes.  The  ground  wave  should  be  free  of  polarization 
changes,  but  this  is  not  the  case  for  the  sky-wave.  The  largest  error  dye 
to  the  sky-wave  occurs  when  the  sferic  originates  at  a  distance  of  about 
300  km,  leading  to  rotation  of  the  magnetic  vector  and  changes  in  the 
apparent  bearing  of  the  source  from  10°  to  30°  (Horner,  1957,  1964).  The  error 
should  diminish  with  decreasing  distance,  however,  and  will  probably  be 
negligible  at  distances  of  interest,  under  200  km,  because  the  ground  wave 
becomes  stronger  at  closer  ranges.  (At  about  300  km,  the  ground-wave  and 
sky-wave  amplitudes  are  comparable.)  One  can  also  "gate"  the  initial  part 
of  the  signal  to  discriminate  against  the  sky-wave  which  is  delayed  in  time 
(W.  L.  Taylor,  private  communication,  1980). 


3.  Range  by  Multiple  Stations 

To  obtain  the  range  of  a  lightning  discharge,  two  crossed-loop  bearing- 
detector  stations  may  be  used  with  triangulation.  This  requires  (a)  that 
the  baseline  of  the  system  be  of  the  same  order  as  the  range  of  the  source 
(lightning  discharge),  and  (b)  that  the  direction  of  the  source  be  different 
from  that  of  the  baseline.  Otherwise  small  errors  in  reading  the  angle  will 
result  in  large  errors  in  the  determination  of  the  range. 

The  remedy  in  the  case  of  ground-based  stations  is  to  use  3  stations 
at  the  vertices  of  an  equilateral  triangle,  with  sides  (baselines)  on  the 
order  of  tens  of  kilometers.  This  arrangement  allows  redundant  verification 
of  the  position  determination  by  each  pair  of  stations.  On  an  airplane  one 
may  be  restricted  to  a  single  station,  and  the  range  would  have  to  be 
determined  by  a  different  means.  Single-station  range  detection  has  proved 
to  be  a  difficult  problem.  A  number  of  solutions  have  been  suggested  (see 
subsection  IV.  B.  on  single-station  range  detectors). 
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4.  Stormscope 


A  modern  commercial  version  of  the  crossed  loop,  the  Stormscope  instru¬ 
ment  presently  available  for  airborne  application  (Ryan,  1977),  uses  a  single 
crossed  loop  for  the  determination  of  the  azimuthal  direction  (bearing)  of 
lightning  discharges  with  respect  to  the  aircraft.  The  operating  frequency 
band  is  centered  on  50  kHz,  that  is,  low  frequency.  Ranges  of  40,  100,  and 
200  nautical  miles  may  be  selected.  Lightning  positions  are  displayed  as 
dots  on  a  CRT  "scope".  In  the  present  version  the  dots  do  not  move  so  as  to 
reflect  the  change  in  position  of  the  aircraft.  Such  dots  cannot  show  the 
spatial  extent  of  the  source.  Also,  the  manufacturer  claims  that  Stormscope 
rejects  signals  from  horizontal  discharges.  Hence  some  data  from  nonverti¬ 
cal  discharges  would  be  lost.  The  manufacturer  claims  a  range  accuracy  of 
+10  percent,  minimal  errors  from  reradiation  by  nearby  metal  surfaces,  and 
invulnerability  to  corona  noise  because  of  special  signal  processing  tech¬ 
niques. 

The  Stormscope  is  of  considerable  interest  because  of  its  compact  geom¬ 
etry  and  relatively  low  cost.  Its  accuracy  and  reliability  are  still  in  the 
process  of  being  assessed. 

The  range  is  indicated  by  the  radial  distance  of  the  dot  from  the  center 
of  the  display.  According  to  the  manufacturer,  the  range  of  the  lightning 
discharge  is  obtained  by  computer  evaluation  of  signal  strength,  time  to  peak 
value,  decay  time,  spectral  content,  and  comparison  of  electric  and  magnetic 
field  amplitudes.  The  details  of  the  physical  concept  of  this  evaluation 
could  not  be  obtained  at  the  present  time.  Therefore,  the  accuracy  of  the 
range  determination  is  still  an  open  question. 


i 


5.  Stormscope  Tests 


In-flight  field  tests  of  Stormscope  were  performed  by  the  Air  Force 
Flight  Dynamics  Laboratory  in  1977  (Seymour  and  Baum,  1979;  Baum  and  Sey¬ 
mour,  1980)  with  funding  by  the  FAA.  (See  also  Rozelle,  1979.)  The 
instrument  was  installed  near  the  leading  edge  of  (and  on  the  underside  of) 
the  right  wing  tip  of  a  T-39B  aircraft.  This  location  was  chosen  to  reduce 
spurious  responses  because  the  instrument  is  very  sensitive  to  electrical 
noise  generated  within  the  aircraft. 

The  tests  were  performed  near  thunderstorms  at  Kennedy  Space  Center, 
in  conjunction  with  on-board  air  weather  radar  (X-band),  ground  weather 
radar,  and  the  ground-based  LDAR  system  for  locating  lightning  discharge 
sources.  (LDAR  is  an  elaborate  time-of-arrival  system  discussed  later.) 

Some  key  results  of  the  Stormscope  tests  are: 

1)  Some  reasonable  correlations  are  obtained  of  Stormscope  dots 
with  weather  radar  precipitation  echoes,  with  some  discrepan¬ 
cies. 

2)  Stormscope  dots  are  much  more  dispersed  both  in  range  and 
azimuth  than  LDAR  dots,  although  roughly  in  the  same  general 
area. 

3)  Isolated  Stormscope  dots  sometimes  appear  in  the  fair- 
weather  region  roughly  180°  opposite  the  main  activity. 

We  believe  these  could  be  due  to  the  180°  sensing 
errors  discussed  above. 

4)  With  respect  to  range,  Stormscope  tends  to  indicate  activity 
more  distant  than  LDAR.  Moreover,  Stormscope  sometimes  shows 
a  "radial  spread"  or  "spoke"  effect,  a  radial  distortion  of 
dots  apparently  due  to  range  errors. 

In  summary  the  authors  (Baum  and  Seymour)  feel  that  despite  the  apparent 
over-spread  of  its  dots  (which  probably  results  from  range  and  bearing  errors), 
Stormscope  provides  a  conservative  warning  of  severe  weather,  i.e.  larger 
avoidance  areas  than  LDAR. 

Spurious  indications  by  Stormscope  (high  "false  alarm  rates")  have  also 
been  observed  in  ground-based  tests  for  mining  applications  (Johnson  et  al ,  1980). 

On  the  other  hand,  it  should  be  noted  that  some  private  pilots  have 
found  Stormscope  helpful  in  avoiding  thunderstorms  (R.  Rozelle,  R.  Collins, 
personal  communications, 1980). 
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6.  Gated  Wide-Band  Direction-Finding  System 

A  recently-developed  crossed-loop  detector  system,  called  the  Lightning 
Location  and  Protection  (LLP)  system,  operates  with  gated  wid'e-band  (GWB) 
electronics  (Krider  et  al ,  1976,  1980).  This  system  is  presently  ground-based 
and  is  designed  to  locate  exclusively  cloud-to-ground  return  strokes,  by  tri¬ 
angulation  using  two  or  more  well -separated  stations.  It  rejects  signals  due 
to  intracloud  discharges,  background  noise,  and  return  strokes  lowering  posi¬ 
tive  charge  to  ground.  For  an  accepted  signal  the  bearing  error  due  to  slant 
in  the  lightning  channel  can  be  reduced,  to  under  1°  according  to  NSSL  tests 
by  M.  W.  Maier  (personal  communication,  1980),  by  gating  on  the  first  few 
microseconds  of  the  signal.  This  corresponds  to  radiation  from  the  lowest 
100  m  of  the  return  stroke  channel.  The  misdirection  is  minimized  because  the 
source  and  the  detector  are  both  essentially  on  the  ground  plane  (see  "Slanted 
Lightning  Channel").  The  LLP  system  uses  1  kHz  to  1  MHz  as  its  operating  fre¬ 
quency  band.  It  is  commercially  available  and  has  recently  seen  widespread 
ground-based  deployment  by  the  USDI  Bureau  of  Land  Management  to  detect  forest 
fires  (Krider  et  al ,  1980).  The  LLP  system  can  be  used  as  close  as  1  km  from 
a  return  stroke  channel  (M.  A.  Oman,  personal  comuni cati on) .  A  similar  sys¬ 
tem  has  been  reported  by  Bent  (1979),  using  however  1  kHz  to  100  MHz. 

With  respect  to  possible  airborne  application,  a  single-station  version 
of  the  LLP  system  probably  can  be  adapted  to  locate  return  strokes  and  display 
the  data  in  a  manner  similar  to  Stormscope.  The  LLP  system  would  as  well  as 
Stormscope  be  subject  to  site  errors.  These  can  be  minimized  by  appropriate 
calibration. 

7.  Kohl's  Multiple  Loops 

A  system  of  multiple  magnetic  loops  has  been  developed  by  Kohl  (1966) 
for  determining  the  direction  of  sources  of  RF  signals  including  sferics. 

It  is  in  essence  a  removal  of  the  two-component  limitation  of  the  crossed 
loop  antenna,  and  extension  to  a  3-component  sensor  by  the  addition  of  a 
third  loop  in  the  third  direction.  Airborne  experiments  have  been  spon¬ 
sored  by  AFGL,  and  the  system  has  been  flown  experimentally  on  a  Canadian 
ASW  patrol  aircraft  (D.A.  Kohl,  private  conmuni cation,  1980).  The  system 
operates  at  500  kHz  and  is  capable  of  detecting  a  magnetic  vector  of  arbi¬ 
trary  orientation,  using  appropriate  electronics.  It  therefore  has  in 
principle  a  (vector  detection)  capability  greater  than  that  of  a  single 
pair  of  crossed  loops.  This  implies  for  example  that  bearing  errors  due  to 
slanted  lightning  channels  can  be  eliminated. 


-TESTS 
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F..  RF  Sferics  Single-Station  Range  Detectors 

It  was  mentioned  earlier  that  single-station  range  detection  is  a 
difficult  problem  area.  Some  proposed  solutions  will  be  discussed  in 
this  section. 


1 .  Stormscope 

In  our  discussion  of  Stormscope  we  referred  to  its  range  indication. 

One  would  infer  from  the  patent  description  (Ryan,  1977)  that  the  range 
is  obtained  from  assuming  that  the  detected  magnetic  field  intensity  H  is 
inversely  proportional  to  the  distance  r.  If  one  is  sufficiently  far  from 
a  radiating  dipole  (so  that  the  induction  term  is  negligible  compared  with 
the  radiation  term  -  see  Eq.  (9))  then  H=K/r,  where  K  is  a  constant  propor¬ 
tional  to  the  time-derivative  of  the  current.  If  K  were  known  and  truly 
constant  (all  lightning  discharges  having  the  same  di/dt),  then  the  range 
would  be  simply  given  by  K/H. 

However,  lightning  discharges  are  highly  variable  and  are  distributed 
with  respect  to  di/dt.*  Data  on  parameters  of  lightning  discharges  are 
difficult  to  obtain.  For  cloud-to-ground  strokes,  data  have  been  obtained 
by  Berger.  Some  of  this  is  tabulated  by  Golde  (1977).  The  data  for  nega¬ 
tive  strokes  may  be  fitted  roughly  by  a  log-normal  distribution  with  stan¬ 
dard  deviation  about  5  dB.  The  corresponding  figure  given  by  Pierce  (1977b) 
is  7  dB  for  VLF  emissions.  This  standard  deviation,  corresponding  to  fac¬ 
tors  of  about  2  and  1/2  with  respect  to  the  median,  would  thus  apply  also 
to  the  range  determinations.  Hence  variability  alone  would  lead  to  both 
underestimates  and  overestimates  in  range  by  factors  of  2. 

An  additional  problem  with  using  magnetic  field  intensities  is  site 
error.  As  shown  in  App.  B,  depending  on  the  location  of  the  instrument 
on  the  airplane,  the  apparent  magnetic  field  intensity  H  can  be  enhanced 
by  as  much  as  a  factor  of  2.  Under  the  r*K/H  assumption,  this  would  lead 
to  underestimates  of  r  by  a  factor  of  2.  The  site  effect  can  also  cause 
diminutions  of  H,  leading  to  overestimates  of  r. 


This  may  contribute  to  the  Stormscope  "spoke"  effect  (see  p.  40). 
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2.  Kohl's  500-kHz  Sferics  Range  Detector 

Fo".  lowing  investigations  by  Nov'inder,  Mai  an,  and  Horner  of  sferics 
radiation  properties,  Kon!  (’.969,  and  references  cited  therein)  notes  that  at 
5Q0  kriz  most  of  the  radiation  is  not  from  the  return  stroke  per  se  but 
rather  from  the  oreakdown  processes  associated  with  the  lightning  (e.g. 
predischarges) .  Using  a  simple  omnidirectional  electric  antenna  (wire), 
the  pulse  spectrum  at  500  kHz  is  found  to  contain  maximum  pulse  amplitudes 
essentially  the  same  for  all  lightning  bursts.  Thus  peak  pulse  measure¬ 
ments  using  spectrum  analyzers  can  form  the  basis  of  a  range  determination 
when  oroperly  calibrated  against  radar  echoes.  The  variation  is  monotonic 
with  range  and  occurs  from  direct-path  propagation  losses  (i.e.,  ground 
waves).  The  variation  depends  on  the  conductivity  of  the  soil. 

Empirically,  one  must  detect  pulses  over  several  minutes  to  have  a 
h i gn  proDability  of  having  detected  the  peak  pulse(s).  That  is,  large 
numbers  of  pulses  ('500)  must  be  detected.  Also,  large  amounts  of  radar 
data  are  required  to  reduce  the  probable  error.  Kohl  (1969)  reports  a 
standard  range  detection  deviation  of  +2.6  Km  over  all  data  ranging  from 
2-274  km.  In  a  recent  private  communication  he  claims  his  range  error  to 
oe  under  10  percent  in  the  range  25-200  miles.  In  addition,  a  rough  fit 
to  his  data  is  given  by  amplitude  proportional  to  l/r'*\ 

A  difficulty  of  the  calibration  with  respect  to  applicability  to  air¬ 
borne  use  arises  because  a  new  calibration  is  initially  required  at  any 
new  location  because  of  the  changed  ground  conductivity  characteristics. 

Regarding  the  existence  of  an  invariant  naximum  pulse  height.  Kohl  says 
that  there  is  (at  500  kHz)  "an  apparent  maximum  limit  of  the  radiation 
energy  at  tnat  frequency.  The  nature  of  the  sources  at  the  instant  that  the 
maximum  occurs  is  ur.  mown  but  it  arises  during  the  complicated  electrical 
breakdown  processes  associated  with  lightning  (the  growth  of  the  leader). 

This  limit  appears  to  be  independent  of  major  variations  in  return  strokes." 


•  % 
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3.  Spectral  Amplitude  Ratio  and  Group  Time  Delay  Difference 

There  is  a  possibility  of  using  spectral  amplitude  ratios  (SAR), 
employing  two  VLF  frequencies  and  sharply-tuned  receivers,  for  sferics 
range  determination  using  propagation  characteristics  of  the  earth- 
ionosphere  wave  guide.  Thus,  for  example,  Schafer  et  al  (1980)  use  two 
frequencies,  5  kHz  and  9  kHz  and  employ  a  propagation  model  to  obtain  a 
linear  relation  relating  the  spectral  amplitude  ratio  SAR  (ratio  of  volt¬ 
age  outputs  at  the  two  frequencies  from  a  whip  antenna)  to  range  R: 

SAR=A+BR,  where  A  and  B  are  constants  given  by  the  model.  They  assume  the 
far  field  and  a  perfectly  conducting  earth,  plus  the  "first  mode"  approxi¬ 
mation.  Hence,  the  range  R  is  obtained  when  SAR  is  measured.  Use  of  SAR 
is  also  suggested  by  Pierce  (1956). 

Another  parameter  is  also  used  based  on  the  high  dispersion  character¬ 
istics  of  the  waveguide  in  the  range  3-10  kHz.  This  is  the  group  delay 
time  difference  (6DD),  which  is  defined  from  the  phases  of  3  equidistant 
frequency  groups,  at  5,  7,  and  9  kHz.  The  GDD  is  formed  from  the  propaga¬ 
tion  times  of  the  groups  at  6  and  8  kHz  (5  to  7  and  7  to  9,  respectively). 
From  the  theory  of  VLF  propagation,  assuming  the  first  mode  only  yields 
GDD=KR'  so  that  R  is  determined  when  the  GDD  is  measured.  According  to 
Pierce  (1977a)  the  GDD  parameter  is  presently  preferred  for  estimating  dis¬ 
tance. 

In  the  near  field  (less  than  300  km  range)  the  interference  of  several 
higher  order  modes  in  the  propagation  of  VLF  waves  gives  rise  to  ambiguities 
in  the  data.  Hence  SAR  and  GDD  data  are  used  in  conjunction.  The  minimum 
range  such  that  the  higher-order  mode  interference  may  be  neglected  is  about 
300  km. 

It  may  be  possible  to  adapt  this  method  to  airborne  use  (H.  Volland, 
private  communication,  1980).  However,  the  method  seems  to  require  large 
amounts  of  statistical  data,  which  is  time-consuming.  Hence  adaptation 
would  require  careful  consideration. 
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Other  possible  range-detection  methods  can  make  use  of  the  dependence 
on  distance  of  electric-field  amplitude  (E)  or  the  ratio  of  magnetic-to- 
electric  field  amplitudes  (H/E).  In  order  to  clarify  the  relationship  of 
field  amplitude  to  source  strength  and  distance,  we  model  the  source  by 
an  oscillating  dipole  and  consider  the  equations  for  the  electric  and  mag¬ 
netic  fields: 

E(r,t)  =  M(t)/4irer3 

+  dM(  t) /dt/4i;ecr2 

+  d2M(  t) /dt2/4-?reC2r  (8) 


and 

H(r,t)  =  udM(t)/dt/4irr2 

+  ud2M(t)/dt2/4ircr  (9) 

Here  we  denote: 

E(r,t),  H(r,t)  =  electric  and  magnetic  field  amplitudes 
M(t)  =  dipole  moment  (charge  x  length) 

r  =  distance  of  receiving  station  from  centerpoint  of  dipole 
c  =  velocity  of  1 ight 
e  =  electric  permittivity  of  air 
u  =  magnetic  permeability  of  air 

These  equations  (which  can  be  found  in  equivalent  form  in  many  textbooks  (e.g. 
Pierce,  1977a))  are  based  on  the  following  assumptions:  The  radiating 
source  is  a  dipole  of  length  L  with  its  axis  vertical.  Its  length  and 
orientation  are  constant  in  time,  and  the  current  is  constant  along  its 
length.  The  distance  r  is  large  compared  with  l.  (Lightning  discharges  do 
not  generally  fulfill  all  of  these  requirements.) 

The  equations  are  valid  for  a  receiving  station  at  the  ground  where 
the  electric  field  has  only  one  component,  perpendicular  to  the  ground 
plane,  while  the  magnetic  field  has  one  (azimuthal)  component  parallel  to 
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the  ground  plane.  Therefore  it  is  possible  to  express  E  and  H  as  scalar 
functions.  For  an  observer  in  space  above  the  ground  plane  (e.g.  on  an 
airplane)  If  and  it  are  three-dimensional  vectors,  and  Eqs.  (8)  and  (9)  must 
be  supplemented  each  by  two  more  equations  for  the  other  two  vector  compo¬ 
nents  . 

3 

In  Eq.  (8)  the  1/r  term  is  designated  as  the  "electrostatic  field" 

2 

term,  while  in  (8)  and  (9)  the  1/r  and  1/r  terms  are  designated  as  the 

3 

"induction  field"  and  "radiation  field"  terms,  respectively.  The  1/r 
term  dominates  in  the  "near  zone,"  while  the  1/r  terms  dominate  in  the 
"far  zone." 

The  values  of  the  scalars  (or  vectors)  E  and  H  depend  not  only  on  range 

r  but  also  on  the  time-derivatives  of  the  dipole  moment,  M,  dM/dt  and 
2  2 

d  M/dt  .  If  the  electric  (or  magnetic)  field  amplitude  is  to  be  used  as  an 
indicator  of  range,  it  must  be  assumed  that  M  and  its  time-derivatives  are 
the  same  -  or  approximately  the  same  -  for  all  lightning  discharges 
detected.  This  is  not  the  case,  however,  and  variability  among  different 
lightning  discharges  can  amount  to  an  order  of  magnitude  or  more.  A  nearby 
weak  discharge  may  thus  produce  a  field  strength  of  magnitude  similar  to 
that  of  a  distant  strong  discharge. 

Therefore  a  dispersion  or  spread  may  be  expected  in  data  representing 
field  amplitude  versus  distance.  The  data  in  Fig.  7,  representing  electric 
field  amplitude  versus  distance  (Kasemir,  1962)  illustrates  such  a  spread. 

The  electric  field  and  distance  are  given  by  the  ordinate  and  aoscissa, 
respectively.  The  field  was  measured  with  a  capacitive-loaded  antenna  with 
a  time  constant  of  6  sec.  The  distance  was  obtained  by  tri angul ation  from 
visual  observations  at  two  stations  separated  by  a  baseline  23.5  km  in 
length. 

We  may  read  this  figure  in  two  ways.  One  way  is  to  determine  the 
spread  in  field  values  at  a  given  distance,  (e.g.  30  V/m  to  230  V/m  at  20  km, 
and  3.2  V/m  to  10  V/m  at  70  km).:'  We  may  on  the  other  hand  determine  the 
spread  in  distance  values  at  a  given  field  value,  e.g.  15  km  to  23  km  at 
200  V/m,  and  66  km  to  90  km  at  4  V/m.  The  first  set  of  data  is  more  nearly 
representative  of  the  "near  zone"  given  by  the  first  term  on  the  right-hand 
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side  of  (8),  and  has  a  range  spread  of  +21%  about  its  central  value  of 
19  km.  The  second  set  of  data  is  more  nearly  representative  of  the  "far 
zone"  given  by  the  third  term  on  the  right-hand  side  of  (8),  and  has  a 
range  spread  of  +19%  about  its  central  value  of  78  km.  It  appears  from 
these  and  other  samples  that  the  relative  error  is  about  +20%  of  the  range, 
independent  of  the  range  itself.  This  gives  a  rough  idea  of  the  error  to 
be  expected  if  the  electric  field  amplitude  is  used  for  range  detection. 

3 

The  solid  curve  in  Fig.  7  corresponds  to  a  1/r  variation.  For  comparison, 

an  inverse-distance  law  (dashed  curve)  is  also  shown,  normalized  to  the 

solid  curve  at  100  km  and  2  V/m.  It  is  clear  from  these  results  that  the 

1/r3  curve  fits  the  data  quite  well  almost  to  100  km.  However,  the  dashed 

curve  indicates  that  a  1/r  variation  better  represents  the  data  beyond 

about  80  km.  It  is  suggested  that  a  fit  to  the  data  may  be  more  effective 

3  2 

as  well  as  physically  more  nearly  correct  using  3  terms  (1/r  ,  1/r  and  1/r) 
corresponding  to  Eq.  (8). 

O 

Range  detection  based  on  the  l/rJ  assumption  has  also  been  investigated 
by  Ruhnke  (1962).  Both  his  results  and  data  of  the  same  type  obtained  earlier 
by  Pierce  (1955)  appear  to  show  a  considerably  greater  spread  than  that  indi¬ 
cated  by  Fig.  7. 


5.  Range  Oeteqtion  by  Flash  Counters 

A  related  problem  arises  in  the  error  determination  in  the  measure¬ 
ment  of  the  lightning  density  (flashes  per  square  kilometer)  by  lightning 
flash  counters  (Horner,  1960).  That  is,  one  records  the  discharges  within 
a  defined  radius  from  the  instrument.  A  wide-band  low-frequency  receiver 
is  used,  responding  to  sferics  electric  fields  having  amplitudes  greater 
than  a  certain  threshold  level.  The  problem  is  then  to  determine  the  dis¬ 
tance  that  corresponds  to  this  threshold  value. 

It  was  discovered  that  the  accuracy  of  the  radius  determination 
was  poor  due  to  the  great  variability  of  the  source  strength.  However,  by 
suitable  statistical  treatment  of  large  amounts  of  data,  Horner  (1960) 
establisf.r'l  an  effective  range  by  comparisons  between  counter  records  and 
observations  by  visual  and  aural  techniques. 


V  • 
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In  addition,  there  is  a  strong  influence  on  the  received  signal  strength 
by  the  differences  in  wave  propagation  characteristics  during  day  and  night. 
This  effect  is  usually  ignored  by  reports  on  testing  of  range  indicators 
based  on  the  amplitude  of  lightning-generated  electromagnetic  pulses.  Tests 
for  airborne  applications  of  these  techniques  should  include  both  day  and 
night  data. 


6.  Range  Detection  by  H/E  Ratio 

To  eliminate  the  influence  of  the  source  strength  represented  by  the 
dipole  moment  and  its  time-derivatives  in  (8)  and  (9)  Ruhnke  (1971a)  sug¬ 
gested  the  use  of  the  ratio  of  the  magnetic  to  electric  field  amplitudes 

H/E.  In  the  near  zone  this  ratio  increases  linearly  with  r  because  of  the 
3  2 

1 / r  dependence  of  E  and  the  1/r  dependence  of  H.  In  the  far  zone  the 
ratio  becomes  constant  since  E  and  H  are  both  proportional  to  1/r.  The 
problem  here  is  that  in  the  near  zone,  where  there  is  a  variation  of  H/E 
with  r,  E  is  proportional  to  M  while  H  is  proportional  to  dM/dt.  Thus, 
for  range  detection  it  must  be  assumed  that  (dM/dt)/M  is  constant  for  all 
discharges.  This  implies  that  M(t)  must  be  exponential  in  time.  As  yet 
there  is  no  evidence  that  this  holds  even  approximately.  However,  it  would 
seem  worthwhile  to  pursue  this  approach.  One  would  have  to  take  into 
account  also  the  day  and  night  differences,  or  assume  that  these  are  elim¬ 
inated  by  the  same  method. 
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C.  Special  Systems 

In  this  section  we  consider  some  examples  of  special  systems  for 
direction-finding  and  location  of  sferics  sources.  Time-of-arri val  (TOA) 
systems,  an  interferometer  system,  and  an  electrograph  system  are  included. 


1.  Time-of-Arrival  (TOA)  Systems 

Two  examples  of  Time-of-Arrival  (TOA)  systems  for  detecting  sferics 
from  discharges  in  electrified  clouds  are 

(a)  LDAR  (Lightning  Detection  and  Ranging)  (Poehler  and 
Lennon,  1978  and  1979) 

(b)  Taylor's  Lightning  Mapping  System  (Taylor,  1978) 

The  two  systems  operate  similarly  in  that  all  the  antennas  in  an  array 
receive  pulse  signals  radiated  from  the  same  discharge  source  located  at  a 
point  in  space  at  a  given  time.  Using  the  differences  in  the  times  of 
arrival  of  the  pulse  at  the  individual  antennas,  an  analysis  may  be  performed 
yielding  the  coordinates  of  the  source  (angles  of  azimuth  and  elevation,  plus 
range  with  a  sufficient  number  of  antennas). 

The  analysis  generally  depends  on  the  fact  that,  for  any  one  pair  of 
antennas  and  a  given  difference  in  TOA  of  the  pulses  to  them,  the  source  must 
be  located  somewhere  on  a  unique  hyperboloid  of  revolution  with  the  two 
antenna  positions  at  the  focii.  For  another  pair  of  antennas  in  the  group 
detecting  the  same  pulse,  there  is  another  hyperboloid  of  revolution.  The 
intersection  of  the  two  hyperboloids  gives  a  curve  in  space  on  which  the 
source  must  lie.  A  third  hyperboloid,  associated  with  a  third  pair  of 
antennas,  intersecting  with  the  latter  curve  serves  to  locate  the  source 
position. 

The  accuracy  of  source  location  in  space  should  be  high  if  the  source 
is  within  a  range  comparable  to  the  effective  baseline  dimension  (i.e., 
average  baseline  dimension  for  a  group  of  several  antennas).  If  the  source 
is  far  away  compared  with  the  effective  baseline  dimension,  the  range  is 
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difficult  if  not  impossible  to  determine,  but  the  direction  (azimuth  and 
elevation  angles)  is  determinable.  In  this  limit  the  hyperboloids  become 
cones,  and  the  intersections  become  straight  lines.  It  should  be  noted  that 
using  two  or  more  small-baseline  groups  of  antennas  with  the  groups  separated 
by  large  distances  increases  the  effective  baseline  dimension.  Thus,  the 
LDAR  system  typically  uses  4  stations  in  a  Y-configuration  with  an  effective 
baseline  of  10-20  km,  so  that  the  range  of  sources  within  10-20  km  of  the 
center  can  be  determined  accurately.  The  Taylor  system,  on  the  other  hand, 
uses  two  groups  of  antennas,  separated  by  about  20  km,  with  baselines  of 
about  14  m  within  each  group. 

Both  of  the  above  T0A  systems  operate  in  the  approximate  frequency 
range  20-80  Wiz.  A  small  amount  of  radiation  from  return  strokes 
(<  10  MHz)  is  detected.  Both  systems  also  depend  on  fast  rise  times  of 
the  signals  to  mark  the  times  of  arrival.  Hence  weak  signals  and  slow 
rises  are  rejected. 

Two  kinds  of  time  resolution  are  involved: 

1)  The  time  resolution  for  a  given  pulse,  which  controls  the 
accuracy  of  the  angle  determination,  and 

2)  The  time  window  for  resolving  individual  pulses. 

The  time  window  for  resolving  individual  pulses  appears  to  be  about  50ns 
for  Taylor's  system,  and  about  lOOvs  for  LDAR.  These  intervals  are  dic¬ 
tated  by  the  systems  used  for  processing  the  signals.  The  T0A  systems  can 
typically  detect  50-100  discharges  accompanying  a  ground  strike,  occurring 
within  milliseconds  of  the  strike.  Reception  rate  (e.g.  25,000/s)  is  inde¬ 
pendent  of  the  time  window  for  receiving  individual  pulses  (W.  L.  Taylor, 
private  communication,  1980). 

With  respect  to  the  time  resolution  for  a  given  pulse,  and  its  angle 
determination,  the  angular  error  (Ae)  is  related  to  the  time  resolution 
(At)  approximately  by 

A6°~  60^  (10) 

where  c  is  the  speed  of  light,  and  d  is  the  baseline  dimension.  One  method 
of  estimating  At  is  to  perform  laboratory  tests  on  differences  in  times  of 
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arrival,  if  the  baseline  is  small  enough;  this  involves  careful  adjustment 
of  band-pass,  time-delays,  and  thresholds.  Another  method  is  to  use  simply 
a  number  of  the  order  of  the  inverse  of  the  frequency  bandwidth  of  the  sys¬ 
tem.  The  use  of  different  definitions  may  explain  differences  in  estimates 
of  accuracy.  In  any  event,  the  above  formula  indicates  that  the  LDAR  system 
with  its  long  baseline  (d)  has  a  potentially  high  location  accuracy.  How¬ 
ever,  difficulties  may  arise  because  of  the  long  baseline.  If  the  time 
between  emissions  of  separate  VHF  pulses  in  the  clouds  is  of  the  order  of  or 
less  than  the  pulse  propagation  time  between  stations  (30- 50us)  the  system 
may  become  confused. 

Some  individual  characteristics  and  airborne  applicability  of  the 
two  TOA  systems  are  described  next. 


(a)  LDAR  System 

The  antennas  of  the  LDAR  system  are  separated  by  tens  of  kilometers 
and  are  all  in  the  same  horizontal  (ground)  plane.  This  means  that  the 
source  x-y  position  determination  is  relatively  precise  when  the  source  is 
(well  isolated  in  time  and)  within  a  baseline  dimension  of  the  array. 
Otherwise,  only  the  azimuth  can  be  precise.  Because  of  the  planar  configu¬ 
ration,  elevation  is  inaccurate  at  low  values  (e.g.,  heights  below  300  m). 

A  hyperbolic  direction-finding  system  with  high  accuracy  and  much  larger 
baselines  was  developed  earlier  by  Lewis  et  al  (1962). 

The  extremely  long  LDAR  baseline  dimension  makes  it  obviously  inapplica¬ 
ble  as  an  airborne  system.  Nevertheless,  this  instrument  has  its  value  for 
testing  airborne  systems  (e.g.,  Stormscope  versus  LDAR  tests),  and  it  can 
also  potentially  be  useful  for  communicating  lightning  activity  data  to 
pilots  from  the  ground. 

(b)  Taylor's  Lightning  Mapping  System 

As  opposed  to  using  arrays  of  antennas  in  the  ground  plane  as  in  LDAR, 
Taylor's  group  consists  of  two  pairs  of  antennas,  one  with  a  horizontal 
baseline,  and  the  other  with  a  vertical  baseline.  The  horizontal -baseline 
pair  is  used  for  azimuth  alone,  and  the  vertical -baseline  pair  is  used  for 
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elevation  alone.  The  baseline  of  this  group  is  of  the  order  of  10  m  (as 
opposed  to  the  10-km  scale  of  the  LDAR  group).  Thus,  its  output  is  azimuth 
and  elevation  (a  direction  line).  Taylor  uses  two  such  groups,  separated 
on  the  order  of  10  km,  plus  triangulation  to  determine  range. 

The  relatively  small  baseline  of  the  Taylor  group  seems  more  applicable 
to  airborne  use  than  the  LDAR.  However,  the  reduction  of  the  baseline  to  the 
order  of  10  m  or  less  (suitable  for  aircraft)  can  result  in  large  angular 
errors.  For  example,  if  we  use  Eq.  (10)  with  d=3.6  m,  and  At  —(60  MHz)”^ 
-1.7x10"^  sec  from  bandwidth  considerations,  we  obtain  A0=85°.  This  result 
suggests  that  high  precision  may  be  difficult  to  obtain  with  TOA-type  air¬ 
borne  systems.  However,  the  prediction  of  errors  is  apparently  not  as 
straightforward  as  implied  in  the  foregoing.  In  this  connection,  Taylor 
(private  communication,  1980)  states  the  following: 

"Bandwidth  considerations  for  obtaining  At  are  totally  irrelevant  for 
Taylor's  system  since  the  difference  in  time  of  arrival  at  each  pair  of 
antennas  is  obtained  from  the  initial  rise  of  the  pulse  waveform.  Labora¬ 
tory  tests  showed  At  could  be  measured  to  <0.5ns  (no  shorter  interval  was 
attempted).  Field  tests  using  the  LDAR  calibration  pulse  transmitter  atop 
the  VAB  (Vehicle  Assembly  Building)  at  KSC  and  using  an  airborne  pulse 
transmitter  on  a  NASA  aircraft  for  Oklahoma  calibrations  indicated  the  azi¬ 
muth  and  elevation  angles  were  determined  to  <0.5°  error,  thus  showing  that 
At  was  obtained  in  an  operational  mode  to  <0.4ns.  For  d=l  m  in  Eq.  (10) 
we  therefore  have  A0=7.2°.  For  d=3.6  m,  very  reasonable  for  aircraft 
installation,  we  have  Ae=2°." 


2.  Interferometer  System 

This  instrument  has  been  used  in  radio-astronomy  to  accurately  locate 
extraterrestrial  sources  of  ratio  emission.  It  has  been  adapted  to  lightning 
location  by  Warwick  and  Hayenga  (Warwick  et  al ,  1979;  Hayenga,  1979).  In  the 
present  version,  a  10-percent-relati ve-bandwidth  receiver  receives  VHF 
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radiation  at  34.3  MHz  emitted  by  breakdown  processes  occurring  at  the  tip 
of  a  lightning  channel.  The  frequency  f  is  one  of  several  adjustable 
parameters,  to  be  defined  later.  The  relative  phases  of  the  signals  arriv¬ 
ing  at  a  pair  of  omnidirectional  antennas  contain  the  desired  information 
regarding  the  direction  of  arrival  of  the  signals.  The  antennas  are 
separated  in  the  current  version  by  a  baseline  d  equal  to  twice  the  wavelength 
x.  (In  the  present  version,  d=17.4m=2X).  The  accuracy  of  determination  of 
the  source  direction  depends  on  the  accuracy  with  which  the  relative  phase 
can  be  determined.  The  determination  is  simplified  by  mixing  the  outputs  of 
the  antennas  with  local  oscillator  signals  offset  in  frequency  by  an  amount 
fQ  (typically  2Q0  kHz)  much  lower  than  f.  The  signals  are  then  multiplied 
with  each  other,  producing  an  interference  pattern  with  a  sinusoidal  modula¬ 
tion  having  the  frequency  f  .  The  phase  of  the  modulation,  which  can  be 
determined  accurately  from  successive  zero-crossing  times  of  the  signal,  is 
directly  related  to  the  relative  phase  of  the  original -frequency  signals 
arriving  at  the  two  antennas. 

The  inherent  high  accuracy  of  source-direction  determination  by  this 
method  is  limited  in  part,  however,  by  the  observation  time  t  (time  of  averag¬ 
ing)  of  a  given  train  of  waves.  The  method  assumes  that  the  train  of  waves 
is  sufficiently  long  to  produce  an  interference  pattern,  and  that  the  radia¬ 
tion  comes  from  a  single  source  of  small  size  during  the  time  of  observation. 
If  the  radiation  comes  from  multiple  sources  the  direction  determination  may 
be  in  error.  Hence  the  received  radiation  is  averaged  over  a  sufficiently 
short  time  interval  to  minimize  the  possibility  of  confusion  with  other 
sources.  In  the  present  design  this  time  interval  x  is  of  the  order  of  l-2us. 
Thus,  pulses  separated  by  2us  or  more  and  their  associated  sources  can  be 
resolved,  and  pulses  lasting  much  longer  than  2ys  can  be  sampled  every 
2ys  to  determine  motion  during  the  pulse. 

From  the  above  phase  shift,  one  infers  the  polar  angle  of  the  source 
with  respect  to  the  baseline  direction.  This  determines  one  angle.  Crossed 
baselines  (two  elements  on  each  line)  give  the  vector  direction  (both  azi¬ 
muth  and  elevation).  Using  two  such  groups  widely  separated  can  yield  the 
source  position  by  triangulation.  The  above  technique  was  designed  (by 
Warwick  ana  Hayenga)  for  detecting  the  steps  of  a  lightning  stepped  leader. 


The  accuracy  of  source-angle  determination  (azimuth  and  elevation) 
depends  on  four  parameters  of  the  system:  wavelength  X,  baseline  d,  band¬ 
width  B,  and  integrating  or  averaging  time  t.  In  degrees,  the  angle  error 
may  be  expressed  (aside  from  a  trigonometric  factor)  as 
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Hence,  for  34.3  MHz  and  d=2X=l7.4  m,  for  B=4  MHz  (roughly  10  percent  of 
the  frequency),  and  for  t-5/B-10’^s ,  we  obtain  the  error: 

A0  -  2.3° 

Another  convenient  design  formula  is  obtained  from  (11)  by  solving  for 
the  frequency  f,  with  B  defined  as  O.lf,  namely, 


f  3/2~ 


90  c 
2 tt  A0  d 


1 


C2t 


(12) 


where  c  is  the  velocity  of  light. 

In  a  possible  application  to  airborne  use,  d  would  decrease.  Let 
d=3.6  m,  as  before.  The  frequency  can  be  increased  to  167  MHz  to  keep 
the  same  ratio  of  x/d.  The  bandwidth  B  would  increase  to  about  16.7  MHz, 
while  the  time  t  can  remain  the  same.  Thus,  Eq.  (11)  would  yield  a  smaller 
A0  because  of  the  larger  value  of  /0T,  yielding  an  even  higher  accuracy  (by 
a  factor  of  about  2). 

The  drawback  to  the  above  set  of  parameters  is  that  the  frequency 
167  MHz  would  contain  considerably  less  lightning  radiation  energy.  The 
lowest  frequency  one  may  use  in  Eq.  (11),  with  d  and  t  fixed,  such  that  a@ 
is  no  larger  than  one  degree,  is  4.6  MHz,  or  X=65  m.  This  frequency  is 
sufficiently  low  to  contain  significant  energy.  However,  there  may  be 
considerable  extraneous  broadcast  noise  in  this  band. 

The  interferometer  and  TOA  systems  both  seem  to  have  considerable  free¬ 
dom  in  the  choices  of  the  parameters.  They  seem  readily  adaptable  to  air¬ 
borne  use.  The  possibility  of  range  determination  in  addition  to  azimuth 
and  elevation  should  be  looked  into. 
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With  two  stations,  e.g.  one  on  each  wing  of  a  large  airplane,  range 
may  be  determined  by  tri angulation.  The  ranging  accuracy  depends  on  the 
angular  accuracy  of  the  individual  stations.  For  an  interferometer  system 
the  possible  angular  accuracy  increases  with  the  (square  roof  of  the)  inte¬ 
gration  time  t.  While  use  of  as-smal 1-as-possible  values  of  t  (e.g.  2us) 
would  allow  finer  details  of  discharges  to  be  detected,  the  resulting  large 
amounts  of  data  would  be  more  appropriate  for  a  ground-based  research  sys¬ 
tem  than  for  an  airborne  lightning  warning  system.  An  interferometer  sys¬ 
tem  with  increased  t  (more  accuracy,  less  detail)  is  relatively  easy  to 
implement  and  more  suitable  to  airborne  application  (C.  0.  Hayenga,  per¬ 
sonal  communication,  1980).  As  an  example,  it  is  suggested  by  J.  W.  Warwick 
(personal  communication,  1980)  that  (with  two  crossed-base-1 ine  stations  on 
a  large  airplane)  using  t=75  ms  at  600  MHz  can  yield  a  relative  range  accu¬ 
racy  of  33%  at  200  km. 

3.  Electrograph  and  Crossed-Adcock  Antennas 

A  lightning  warning  ( "El ectrograph")  system  operating  at  high  frequen¬ 
cies,  in  the  900  MHz  region,  has  been  proposed  by  E.  A.  Lewis  and  his  co¬ 
workers.  A  directional  parabolic  antenna  scans  both  azimuth  and  elevation 
to  detect  radio  noise  pulses  originating  in  small-scale  electrical  dis¬ 
charges  ( pre-di scnarges)  preceding  overt  lightning  flashes  (Harvey  and 
Lewis,  1972).  In  this  sense  the  system  is  similar  to  the  lightning-mapping 
T0A  and  interferometer  systems. 

A  sferics  system  for  lightning  bearing  detection  in  the  broadcast  band 
of  frequencies  (0.5  to  1.5  MHz)  was  developed  by  Stergis  and  Doyle  (1959). 
Two  Adcock  antennas  placed  at  90°  to  one  another  were  selected  over  crossed 
loops  because  of  the  Adcock's  lower  sensitivity  to  polarization  errors  (due 
to  down-coming  sky-waves).  Tests  of  crossed  loops  showed  that  20°  to  30° 
bearing  errors  were  not  uncommon.  The  accuracy  claimed  is  about  5°,  with  a 
potential  reduction  to  2°. 

With  respect  to  possible  airborne  applications,  the  Electrograph  is 
being  designed  for  this  purpose  (R.  B.  Harvey,  personal  communication). 

It  is  uncertain  whether  the  Adcock  system  can  be  made  sufficiently 
compact.  Also,  the  Adcock  system  may  be  as  susceptible  as  crossed 
loops  to  aircraft  site  errors.  These  questions  need  investigation. 


V.  OPTICAL  BEARING  DETECTS 

Flash  location  by  using  the  optical  signals  from  lightning  has  some 
similarities  to  the  sferics  methods  discussed  in  Sec.  IV.  For  example, 
Kidder  (1973)  has  described  a  South  African  system  employing  several 
cameras  to  give  bearings,  and  subsequent  triangulation  to  locate  the  dis¬ 
charge  (Pierce,  1977b).  Silicon  photodiode  systems  have  been  used  on  a 
satellite  for  detection  of  lightning  from  space  (Edgar,  1979;  Edgar  and 
Turman,  1980;  Turman  and  Edgar,  1980). 

An  optical  system  for  detection  and  recording  of  lightning  on  the 
ground,  from  aircraft  and  from  space  (the  Orbiter)  is  under  development 
by  Vonnegut  and  his  co-workers  (Griffiths  and  Vonnegut,  1975;  Vonnegut 
and  Passarelli,  1978;  Vonnegut  et  al,  1980).  Their  system  consists  of  a 
photocell  and  a  super-8  sound  motion  picture  camera,  designed  to  be  hand¬ 
held. 

Both  the  Vonnegut  and  the  Edgar-Turman  systems  can  detect  lightning 
in  clouds  illuminated  by  direct  sunlight  even  when  the  lightning  cannot 
be  seen  by  the  human  eye.  Daylight  interference  is  also  avoidable  to  some 
extent  by  use  of  a  filter  for  the  H-alpha  line  which  is  much  stronger  in 
lightning  than  in  daylight. 

Thus  far  the  optical  systems  available  are  suitable  for  picking  up 
signals  from  the  direction  pointed  to  by  the  operator;  they  are  not  auto¬ 
matic  direction-finders  although  various  concepts  for  automatic  systems 
are  presently  under  consideration  (B.  C.  Edgar,  personal  communication, 
1980).  Some  comparisons  of  bearing  detection  of  lightning  between  an 
optical  sensor  and  a  Stormscope  sferics  detector  are  presently  being 
planned  in  terms  of  a  ground-based  study  (0.  H.  Vaughan,  personal  communi¬ 
cation)  . 

With  respect  to  range  detection,  a  limitation  of  an  optical  system 
is  the  line  of  sight.  Hence  such  a  system  would  probably  involve  short- 
range  triangulation. 

With  respect  to  bearing  errors,  in  principle  there  should  be  none  for 
vertical  lightning  channels.  However,  there  is  generally  significant  hori¬ 
zontal  branching,  and  many  lightning  channels  are  approximately  horizontal. 
These  can  represent  broad  sources,  leading  to  large  bearing  errors.  Another 
source-broadening  effect  occurs  even  if  the  channel  is  vertical.  Namely,  the 


light  emitted  by  the  lightning  can  undergo  nearly  isotropic  scattering 
with  little  attenuation  within  the  same  cloud  or  by  nearby  clouds,  so  that  the 
light  source  as  seen  by  the  detector  can  be  broad  in  extent.  Thus,  large 
regions  of  a  cloud  bank  can  be  diffusely  illuminated  and  subtend  a  large 
angle  at  the  detector  because  of  the  light  scattering  (the  whole  cloud 
bank  "lights  up").  For  example,  a  10-km-wide  cloud,  brightly  illuminated 
from  lightning  within  and  at  a  distance  of  50  km  constitutes  a  source  sub¬ 
tending  an  angle  11.4°  at  the  detector.  Thus  the  bearing  error  can  be 
+5.7°. 

In  addition,  the  attenuation  may  become  serious  if  the  line  of  sight 
passes  through  a  large  thickness  of  cloud  material,  e.g.  parallel  to  a 
frontal  system. 

These  questions  need  to  be  considered  in  designing  an  optical  light¬ 
ning  warning  system. 

Some  possible  constructions  of  optical  detectors  are  suggested  by  the 
sketches  in  Fig.  8.  The  sketches  labelled  (a)  in  the  figure  use  the 
dependence  of  photocell  output  on  the  angle  of  incidence  of  the  incoming 
light  (e.g.  a  cosine  law).  The  detectors  on  the  4  side-surfaces  of  the 
cube  or  on  the  8  surfaces  of  the  cross  (E.  P.  Krider,  personal  communica¬ 
tion)  will  have  different  outputs  in  general  dependent  on  the  angle  of  the 
incident  illumination.  These  outputs  can  be  compared  and  analyzed  to  infer 
the  correct  direction. 

The  sketch  labelled  (b)  is  a  rotating  slit  and  mirror,  using  a  single 
photocell.  The  sketches  labelled  (c)  use  a  "fish-eye"  lens  or  a  parabolic 
mirror  to  direct  incoming  light  to  one  of  a  large  number  of  photodiodes 
arranged  in  a  circle. 

The  dimensions  of  the  above  systems  are  very  small,  only  of  the  order 
of  an  inch. 


FIG.  8.  Optical  Detectors 
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VI.  METEOROLOGICAL  INSTRUMENTATION 

Two  airborne  meteorological  instruments  of  interest  in  lightning  warn¬ 
ing  systems  are  the  familiar  airborne  weather  radar  and  a  temperature 
sensor. 

We  consider  first  the  temperature  sensor  which  is  of  interest  because 
the  likelihood  of  encountering  high  electric  fields  and  of  triggering  light¬ 
ning  by  an  aircraft  is  greater  in  some  temperature  regions  than  others  (see 
Sec.  II  on  triggered  lightning  and  statistical  evaluations).  Typical  air¬ 
borne  instruments  for  measuring  temperature  are  a  platinum  wire  and  a 
reverse-flow  temperature  probe  ( Musi  1  and  Prodan,  1980). 


Radar  and  Masking  Effect 


The  airborne  weather  radar  operates  by  reflecting  (scattering)  micro- 
waves  from  precipitation-size  cloud  particles.  It  can  detect  precipitation 
areas  over  a  range  of  perhaps  200  km.  Although  this  system  is  widely 
used,  care  is  required  in  interpreting  the  observed  echoes.  For  example, 
reflection  from  a  nearby  heavy  rain  shower  in  front  of  a  more  distant  severe 
storm  can  mask  the  severe  storm  behind,  and  lead  to  an  erroneous  inter¬ 
pretation  that  the  storm  is  weak.  This  attenuating  factor  depends  on  fre¬ 
quency.  Therefore,  the  design  of  the  system  requires  care.  There  is  less 
attenuation  with  longer  wavelengths  but  then  the  antenna  must  be  larger  to 
keep  the  beam  narrow.  The  standard  wavelengths  used  are  as  follows 
(Battan,  1973) : 


Wavelength  (cm) 
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The  S  and  L  band  radars  are  more  penetrating  than  the  C,  X  or  K. 

However,  they  are  generally  too  large  for  airborne  use  and  are  therefore 
ground-based.  The  S  and  L  bands  would  not  be  subject  to  masking,  but 
would  only  see  large  drops  and  hail,  and  would  therefore  miss  light  showers, 
particularly  at  a  distance.  The  C-band  is  relatively  old  and  has  seen  air¬ 
borne  use,  but  is  being  replaced  by  X-band  units  which  are  lighter  and  less 
expensive.  Although  X-band  radar  will  see  light  precipitation  and  cloud 
particles  down  to  about  100  micrometer  in  diameter,  it  is  subject  to  the 
masking  problem.  Hence,  an  airborne  radar  system  is  a  compromise  between 
range  capability  and  reflectivity  capability. 

The  power  intercepted  by  a  radar  antenna  due  to  back-scattered  radiation 
by  particles  of  linear  dimension  small  compared  with  the  wavelength  is  pro-  I 

portional  to  the  2nd  power  of  the  inverse  wavelength  (Battan,  1973).  With 
respect  to  the  choice  of  radar  wavelengths,  we  quote  Mason  (1971):  "How¬ 
ever,  the  increased  absorption  and  attenuation  of  the  incident  and  scattered 
energy  at  very  short  wavelengths  may  more  than  compensate  for  the  higher 
reflectivity  and  seriously  restrict  the  maximum  range,  and  again  a  compro¬ 
mise  has  to  be  effected.  The  radar  sets  in  general  use  (ground-based)  work 
on  either  10-cm  or  3-cm  wavelengths,  10  cm  being  essential  if  it  is  desired 
to  penetrate  heavy  rain.  Several  people  have  suggested  that  a  wavelength  of 
5.6  cm  would  be  the  most  suitable  for  the  study  of  precipitation,  especially 
snow.  A  radar  working  on  5.6  cm  would  be  roughly  10  dB  more  sensitive  than 
an  otherwise  similar  equioment  working  on  10  cm;  but  for  a  path  of  100  miles 
through  rain  of  intensity  10  mm/h,  the  greater  attenuation  on  5.6  cm  would 
just  offset  this  gain." 

Some  modern  airborne  radar  units  are  digitally  compensated  for  distance, 
and  nave  3  recei ved-power  threshold  settings  which  are  adjusted  for  nominal 
rainfall  rates  of  weak,  moderate  and  severe  (C.  B.  Kitchens,  personal  commu¬ 
nication,  1980). 

Another  type  of  radar  is  doppler  radar  that  detects  the  velocity  of 
motion  of  the  scatterers.  It  potentially  would  be  useful  for  detecting  tur¬ 
bulence  and  strong  updrafts  and  downdrafts.  However,  this  is  not  contem¬ 
plated  for  airborne  use  because  it  would  be  prohibitively  expensive  if  it 
were  intended  solely  for  severe  weather  avoidance. 


1 


Radar  echo  characteristics  have  been  related  to  lightning  by  Fitz¬ 
gerald  (1978).  Cloud- to-ground  lightning  is  most  likely  to  occur  near  the 
edge  of  the  radar  precipitation  cores  and  in  heavy  rain  sheets.  In  cold 
clouds  "it  is  believed  that  dBZ  values  of  25  or  30  represent  probable 
threshold  values  for  initial  cloud-to-ground  lightning  activity....  In 
the  warm  cloud  case,  the  scanty  evidence  available  suggests  35  to  40  dBZ 
as  a  more  likely  radar  threshold." 

Thick  anvils  (overhangs  of  thickness  greater  than  5000  feet)  are 
likely  to  produce  occasional  but  intense  cloud- to-ground  strokes,  although 
they  may  extend  a  number  of  miles  from  the  main  precipitation  cores.  In 
some  cases  one  or  two  isolated  intense  strokes  may  occur  from  the  anvil 
when  the  storm  appears  to  be  dissipating.  These  anvils  are  difficult  to 
detect  with  radar,  "having  only  weak  or  no  radar  return  depending  on  the 
individual  radar  system  characteristics."  They  probably  require  X  or  C 
band  units  (at  close  range)  for  detection. 

Fitzgerald  suggests  that  data  on  time  history  of  echo  intensity  and 
top  height  of  moving  storms  are  particularly  significant  for  radar  light¬ 
ning  warning,  and  “that  it  will  be  easier  to  decide  if  new  cellular 
activity  is  likely  to  produce  lightning  than  it  will  be  to  determine  if 
additional  lightning  will  occur  from  the  older  dissipating  regions  of  a 
storm....  Nearly  all  convective  tops  reaching  -40°C  and  50  dBZ  echoes 
will  produce  some  cloud-ground  lightning." 
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VII.  SUMMARY  AND  RECOMMENDATIONS 

We  have  considered  a  number  of  instruments  for  possible  airborne 
lightning  warning  systems.  All  have  advantages  and  disadvantages.  Some 
are  well  known  and  others  poorly  known. 

In  Table  1  we  assign  some  tentative  ratings  to  the  various  instru¬ 
ments  considered.  They  are  rated  according  to  the  following  categories 
(acronym  column  headings): 

AVAI  =  availability 

SINT  =  simplicity  of  interpretation  of  data 

SCON  =  simplicity  of  construction  (essentially  equivalent 
to  cost) 

SMAI  =  simplicity  of  maintenance 
SCAL  =  simplicity  of  calibration 
EXTT  =  extent  tested  (especially  on  aircraft) 

PACC  =  potential  accuracy 

VERS  =  versatility  (possible  use  for  both  near  and  distant 
warning,  or  for  both  range  and  bearing) 

INSE  =  insusceptibility  to  errors  (e.g.  site  errors)  or  to 

disturbances  associated  with  the  environment  (e.g.  rain) 

The  judgments  are  given  in  terms  of  letters: 

L  =  low 
M  =  medium 
H  =  high 
HH  =  very  high 
U  =  unknown 

Asterisk  (*)  =  in  the  present  state  of  development  -  can  be 
improved. 

The  instruments  are  segregated  into  5  groups,  arranged  vertically 
in  the  table: 

1.  Field  detectors 

2.  Sferics  bearing  detectors 

3.  Sferics  range  detectors 

4.  Optical  bearing  detectors  (one-member  group) 

5.  Airborne  radar  (one-member  group) 
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A.  Individual  Assessments 


(1)  Among  the  field  detectors  we  make  the  following  comments. 

Field  mill :  High  in  all  categories  except  that  they  must  be  constructed, 
and  construction  is  not  simple. 

Radioactive  p-obe:  Similar  to  field  mill,  but  construction  is  much 
simpler.  However,  maintenance  is  not  simple  and  some  errors  are 
possible. 

Corona  point:  Extremely  simple  to  construct  and  maintain,  and  not  sus¬ 
ceptible  to  errors.  Interpretation  may  not  be  simple. 

(2)  Among  the  sferics  bearing  detectors  we  make  the  following  comments. 

Crossed  loops:  Available  and  relatively  simple  to  construct  and  maintain. 
Susceptible  to  site  errors,  requiring  careful  calibration.  Narrowband 
system  tested  with  inconsistent  results.  Potential  accuracy  lower 
than  that  of  gated  wi de-band  system  for  detecting  return  strokes. 

Kohl  multiple  loops:  Must  be  constructed.  Not  extensively  tested. 
Otherwise  similar  to  crossed  loop  except  that  potential  accuracy 
may  be  higher  in  detecting  nonvertical  lightning. 

TOA  and  interferometer:  Must  be  constructed.  Mostly  unknown  properties 
except  that  potential  accuracy  may  be  high. 

(3)  Among  the  sferics  range  detectors,  none  have  as  yet  high  accuracy. 

We  make  the  following  comments. 

Crossed  loops:  Comments  similar  to  those  on  crossed-loop  bearing 

detector.  In  addition,  low  potential  accuracy  and  susceptibility 
to  errors  due  to  lightning  variability  may  be  difficult  to  improve. 

Kohl  500  kHz:  May  be  simple  to  construct,  and  promising,  but  largely 
unknown  characteristics. 

SAR-GDD:  Available,  simple  to  maintain,  and  tested  to  some  extent. 
Otherwise  not  known. 

E-ampl itude:  Available,  versatile,  and  simple  to  interpret,  construct 
and  maintain.  Extensively  tested,  but  usually  the  time  difference 
between  lightning  and  thunder  has  been  used  to  determine  distance, 
and  this  is  not  accurate. 
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H/E  ratio:  Must  be  constructed,  but  simple  to  construct  and  maintain. 
Accuracy  and  error  susceptibility  unknown. 

(4)  Regarding  optical  detectors  and  airborne  radar,  we  make  the  follow¬ 
ing  comments. 

Optical  bearing:  Simple  to  construct,  maintain,  calibrate  and  interpret. 
Must  be  constructed.  Not  versatile.  Subject  to  reflection  (scatter¬ 
ing)  errors  due  to  surrounding  clouds. 

Airborne  radar:  Available  and  extensively  tested.  May  be  difficult  to 
interpret  because  of  possible  masking  effects.  Cannot  positively 
identify  regions  of  electrical  activity  but  valuable  as  an  adjunct  to 
any  of  the  above  systems. 


B .  Recommendation s 


We  believe  that  many  of  these  instruments  have  virtues  worthy  of  con¬ 
sideration  for  airborne  warning  and  avoidance  systems.  However,  they 
require  careful  testing  and  evaluation.  Many  possible  combinations  of 
different  types  of  instrumentation  may  be  effective,  for  example,  field 
and  optical  detectors  in  addition  to  radar.  Also  versatility  of  field 
detectors  with  respect  to  both  near  and  distant  warnings  should  be  con¬ 
sidered,  as  well  as  range  detection  with  multiple  bearing  detectors. 

Analyses  performed  in  this  report  should  be  helpful  in  recognizing 
errors  and  in  suggesting  improvements. 


V  * 
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Table  1.  Characteristics  of  Instruments 


AVAI 

SINT 

SCON 

SMAI 

SCAL 

EXTT 

PACC 

VERS 

INSE 

1.  Field 

Field  Mill 

L 

M* 

L 

H 

H 

H 

H 

H 

H 

Radioact.  Probe 

L 

M* 

H 

L* 

H 

H 

H 

H 

M 

Corona  Point 

L 

L* 

HH 

HH 

M* 

M 

U 

H 

HH 

2.  Sf.  Bearing 

Crossed  Loops 

H 

M* 

H 

H 

L* 

H 

L 

M 

L* 

Kohl  Multiple 
Loop 

L 

M* 

H 

H 

M* 

L 

H 

M 

L* 

TOA 

L 

L 

U 

U 

U 

U 

H 

M 

U 

Interferometer 

L 

L 

U 

u 

U 

U 

H 

M 

U 

3.  Sf.  Range 

Crossed  Loops 

H 

M 

H 

H 

L 

H 

L 

M 

L 

Kohl  500  kHz 

L 

U 

H 

M 

L 

M 

U 

U 

U 

SAR-GDD 

H 

M 

M 

H 

L 

H 

M 

U 

U 

E-amplitude 

(1/r3) 

H1 

H 

H 

H 

M 

M 

L 

H2 

M 

H/E  ratio 

L 

M3 

H 

H 

L 

L 

U 

U 

U 

4.  Optical 

Bearing 

L 

H 

H 

H 

H 

M 

M 

L 

L 

5.  Airborne 

Radar 

H 

M4 

L 

M 

M 

H 

U 

M 

U 

1  Lightning  counter  available 
extended  to  1/r  variation. 

in  near- 

•range 

form,  not 

far 

range. 

Can  be 

2 

Slow  antenna. 


3 

Near-range  only. 

4 

May  be  difficult  to  interpret  when  masking  occurs. 

★ 

In  the  present  state  of  development  -  can  be  improved. 

Notes:  Electrograph  and  crossed-Adcock  antennas  not  included:  character¬ 
istics  uncertain.  LDAR  not  suitable. 

Crossed-loop  assessments  here  refer  to  narrow-band  systems.  A 
gated  wi de-band  system  applied  to  return  strokes  should  be  more 
accurate  for  bearing  detection. 


APPENDIX  A 


Polarization  Errors  in  Lightning  Direction  Finding  by 
Magnetic  Crossed  Loops:  "Misdirection"  Due  to 
Nonvertical  Orientation  of  the  Lightning  Channel 

In  this  section  we  consider  the  error  in  magnetic  direction  finding 
due  to  the  effect  of  nonvertical  orientation  of  a  radiating  lightning 
channel,  modeled  by  a  radiating  tilted  dipole.  The  error  in  direction 
finding  derived  analytically  here  (to  our  knowledge  not  published  previ¬ 
ously)  is  that  due  to  the  tilt  of  the  source  dipole,  where  both  source  and 
observer  are  above  a  conducting  ground  plane.  The  magnetic  crossed-loop 
detector  is  considered  to  be  otherwise  perfect  (no  perturbations  due  to 
the  aircraft,  etc.).  The  "misdirection”  due  to  induced  skin  currents  of 
the  aircraft  (siting  errors)  is  treated  in  Appendix  B.* 

Referring  to  Fig.  A-l,consider  the  detector  to  be  at  altitude  h  above 

the  ground  plane.  (The  loop  is,  of  course,  mounted  on  an  aircraft,  but 

the  perturbing  influence  of  the  aircraft  is  ignored  for  the  present.)  The 

lightning  radiation  source  is  represented  by  a  dipole  of  vector  dipole 

moment  P,  located  at  altitude  H  above  the  ground  plane.  Define  a  Cartesian 

coordinate  system,  wherein  the  ground  plane  is  represented  by  the  x-y  plane, 

and  the  dipole  is  centered  on  the  z-axis  at  a  vertical  height  z=H  above  the 

ground  plane.  The  detector  is  located  in  the  x-z  plane,  at  a  height  z=h 

above  the  ground  plane,  and  at  horizontal  distance  x=D  from  the  z-axis  (for  y=0). 

Thus,  the  coordinates  of  the  source  and  detector  in  this  system  are 

(0,  0,  H)  and  (D,  0,  h),  respectively.  Let  R-j  denote  the  distance 

from  the  dipole  to  the  detector,  and  let  R£  denote  similarly  the 

distance  from  the  dipole  image  (at  depth  H  below  the  ground  plane)  to  the 

same  detector.  Then  R-|  and  R£  are  given  by  R^=D^+(H-h)^  and  R^D^H+h)^. 

Let  the  dipole  be  arbitrarily  oriented,  with  components  P  ,  P  ,  and  P  . 

a  y  z 

We  can  then  compute  the  magnetic  field  as  the  superposition  of  contributions 

*  ~  a  A 

from  the  3  separate  vector  dipoles  of  moments  P  i,  P  j,  and  P  k,  where  i, 

a  a  x  y  z 

j,  and  k  denote  unit  vectors  in  the  x,  y  and  z  directions. 

According  to  Hertz'  solution  of  Maxwell's  equations,  the  electric  and 
magnetic  fields  due  to  a  radiating  dipole  may  be  expressed  in  terms  of  a 
dipole  moment  function  (e.g.,  Sommerfeld,  1952;  Abraham  and  Becker,  1937). 

* 

The  term  "misdirection"  appears  in  the  text  by  Sommerfeld  (1949). 


FIG.A-l.  NONVERTICAL  OIPOLE  ANO  IMAGE 


L» 
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If  we  denote  this  time-dependent  function  by  P(t-R/c),  the  vector  potential 
A  may  be  written  as 


A  -  \  P  (t-R/c) 


0) 


where  the  argument  denotes  the  retarded  time  at  distance  R  from  the  dipole, 

and  the  dot  denotes  differentiation  with  respect  to  the  time.  Here, 

the  units  adopted  are  Gaussian  and  t  is  expressed  in  emu,  i.e.,  abamp-cm, 

with  R  in  cm.  (In  MKS  units,  an  additional  factor  u  /4u  appears  on  the 

0  -7 

right  side,  where  yQ  is  the  magnetic  permeability  of  free  space,  4irxl0 
henry /meter;  then  P  is  expressed  in  ampere-meters, wi th  R  in  meters.) 

Now  the  magnetic  field  intensity  B  is  given  by 

B  =  curl  A(R)  =  (R/R)  x  (dA/dR)  (2) 

since  A  is  a  function  of  R  only.  Differentiation  yields 


-+  -*■  -+ 
dA  _  P  _P_ 
dR  "  '  r2  *  Rc 


(3) 


where  the  second  term  results  from  the  dependence  of  the  time-like  argu 
ment  on  R,  and  we  obtain 


B  =  M(R)  x  R 


(4) 


so  that  a  =  1/R3  and  3  =  1/cR2. 


In  terms  of  our  cartesian  coordinates,  let 


where  we  identify  the  unsubscripted  variables  with  the  original  dipole. 

It  follows  from  (4)  and  (6)  that,  at  the  observer's  position,  the  magneti 
field  of  the  original  dipole  is  given  by 

^  =  t-(H-h)Mly]  i  +  [(H-h)Mlx  +  D  Ml2]  j  +  [-D  Mly]  k  (8) 

while  the  magnetic  field  due  to  the  image  is  given  by 

B2  =  E(H+h)M2y]  i  +  [-(H+h)M2x  +  D  M2z]  j  +  [-0  M2y]  k  (9) 

Hence,  the  sum  of  the  two  fields  is: 
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-+  -> 

B  =  B-,  +  B2 

•  C(M2y  -  VH  +  (M2y  +  Mly)h]  * 

+  C-(M2x-Mlx)H-(M2x  +  Mlx)h  +  (M22+Mu)D]j 

+  C'(M2y  +  VD]  *  0°) 

For  a  vertical  dipole,  all  components  of  and  M2  vanish  except  for 

M-jz  and  M2z>  Then  B  has  only  a  y  component,  and  there  is  no  misdirection. 

When  B  has  an  x-component  also,  the  misdirection  angle  is  given  by  the  arc 

tangent  of  -B  /B  (Fig.  A-lb) .  Consider  the  ratio  -B  /B  ,  obtained  from  (10) 
a  y  x  y 


Using  (5)  and  (7),  we  may  rewrite  (11)  as: 


Bx  pyH  +  qyh 

B  pH  +  q  h  +  p  0 
y  a  x  z. 

where 

px  =  (ct2  +  al)Px  +  (P2  +  B1)Px 

py  =  (a2  +  al)Py  +  (e2  +  el)Py 

Pz  =  (<*2  +  ai)p2  +  (^2  +  &])p2 

qx  =  (a2  -  °l>Px  +  (02  "  el)Px 

qy  =  (“2  "  ai ) Py  +  (e2  '  6l^Py 


(12) 


(13) 


•  • •  •  •*  •  •• 

and  where  P  ,  P  ,  P  ,  P  ,  P  P  refer  to  the  original  dipole  components, 

x  x  y  y  z  9  i 

Divide  through  numerator  and  denominator  by  pzD,  and  note  that  we  may 
wri te : 


\  • 
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P 

x 


P 


y 


p 


(14) 


where  P  ,  P  ,  P  are  the  original  dipole  moment  vector  components,  and 

a  y  £ 

where  P  and  P  are  the  time-derivatives  of  the  vector  magnitude  P. 

Then  Px/Pz,  Py/Pz>  qx/Pz,  and  Py/Pz  become,  using  (13)  and  (14): 


P 

P 


x_ 

z 


P  P 
=  JL 
P_  P_ 


(15) 


where 


( x,  -  a, )  P  +  (  3?  -  3, )  P 

Q  =  — = - 1— ; - - - !— - 

( ot2  +  a-j )  P  +  (  S>2  +  )  'p 

(1/R23  -  I/R^JP  +  (1/R22  -  1/R12)P/c 
(1/R23  +  1  /R-j 3)  P  +  (1/R22  +  1/R12)P/c 


Finally,  the  ratio  -Bx/By  may  be  written: 


P 


1  +  Gpi 
z 


(16) 


(17) 


, _ - -■ 
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where 


_  H  +  Qh 
=  [T^ 


(18) 


If  the  detector  is  on  the  ground  plane  (h=0),  the  problem  simplifies 
greatly,  and  we  have  a-pc^  and  B-|=B2>  so  that  Q=0.  Then  G  reduces  simply 
to  G=H/D,  and  we  have 


H 

D  P 

,  +  H 
1  +  U  PJ 


(19) 


which  is  identical  to  the  result  given  by  Uman  et  al  (1980  -  see  plots),  re¬ 
expressed  in  our  terms,  for  a  detector  on  the  ground  plane.  An  equivalent 
result  was  also  given  earlier  for  this  case  by  Kalakowsky  and  Lewis  (1966). 

If  the  ground  plane  is  completely  nonconducting  (no  image),  then 
the  case  of  the  detector  at  altitude  h  above  ground  also  simplifies 
greatly.  In  this  case  cu  and  62  vanish,  and  Q=-l .  Then 


B 


(H-h)  y 

0  pz 


y  i  +  (H~h)  _x 
'  D  P 


(20) 


That  is,  H  is  replaced  in  (19)  by  the  "net  altitude,"  H-h. 

A  geometric  interpretation  of  (19)  (or  (20))  is  obtained  by  rewriting 
(19)  in  the  form 


(0  +  H  p-^)  8x  +  H  p 


Py 

H  B. 


(21) 
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This  expression  implies  that  the  "misdirection,"  in  the  plane  of  the 

observer  and  perpendicular  to  B  at  the  observer,  passes  through  the  point 

whose  x,y,z  coordinates  are  -(PV/P,)H,  -(PU/P_)H,  and  zero.  This  point 

x  z  y  z 

(indicated  by  the  letter  I  in  Fig.  A-la)  is  also  the  intercept  of  the 
extended  line  of  the  dipole  with  the  z=0  plane.  This  simple  geometric 
interpretation  was  pointed  out  by  Kalakowsky  and  Lewis  (1966). 

In  the  case  of  harmonic  radiation,  of  angular  frequency  <a,  we  may 

•  *  * 

replace  P  by  ju>P  in  (16),  so  that  Q  is  given  by  the  complex  quantity, 
with  k=u>/c: 


(1/R23  -  1/R-j3)  +  j(k/R22  -  k/R^) 
(1/R23  +  1/R-j 3)  +  j(k/R22  +  k/R^) 

(x3  -  x3)  +  jk(x2  -  x2) 

(x3  +  x3)  +  jk(x2  +  x2) 
a  +  jb 
a '  +  jb ' 

where  x  =  1/R  (x^ 

3  3 

a  =  x^  -  x2 

a'=  +  Xj 


j(ab'  -  ba') 
(a1)2  +  (b1)2 


so  that 

(aa'  +  bb ' ) 

Q  =  -  ^  2 

(a1 )2+(b' )2 


=  1/Rj  ,  x2  =  1/R2) 
b  =  k(x2  -  x2) 
b'=  k(x2  +  x2) 


(22) 


(23) 


(24) 


*  Re  Q  +  j  Im  Q 
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APPENDIX  B 

Polarization  Errors  in  Lightning  Direction  Finding  by 
Magnetic  Crossed  loops:  "Misdirection"  Due  to 
Wave  Scattering  by  Nearby  Metal  Surfaces  (Site  Errors) 

In  this  section  we  consider  the  bearing  error  in  magnetic  direction 
finding  due  to  scattering  of  the  incident  electromagnetic  wave  by  con¬ 
ducting  surfaces  near  the  detector.  The  incident  magnetic  vector  is 
assumed  here  to  be  in  the  horizontal  plane  as  produced  by  a  vertical 
current  source.  (Appendix  A  treats  the  misdirection  due  to  a  nonvertical 
current  source.)  This  incident  magnetic  field  induces  eddy  currents  in 
the  nearby  surfaces.  These  currents  in  turn  produce  a  secondary  magnetic 
field  that  results  in  a  total  magnetic  vector  of  generally  different  (dis¬ 
torted)  orientation.  T^e  crossed-loop  detector  is  "misdirected"  when  it 
senses  the  distorted  orientation  and  indicates  an  erroneous  bearing. 

(The  term  "misdirection"  is  used  in  the  text  by  Sotmierfeld,  1949). 

The  following  error  analysis  aoes  not  appear  to  have  been  published 
previously. 

Evaluation  of  the  scattering  of  incident  electromagnetic  waves  is  in 
general  a  difficult  problem,  even  by  perfectly  conducting  bodies  of  sim¬ 
ple  shape,  e.g.  spheres  (Stratton,  1941).  In  principle,  matrix  equations 
of  infinite  order  are  to  be  solved  for  the  expansion  coefficients  (expan¬ 
sion  in  vector  spherical  wave  functions).  The  difficulties  are  greater 
for  spheroids  although  the  vector  wave  equation  is  separable  and  solutions 
are  obtainable,  in  spheroidal  coordinates  ( e.g .,  Taylor,  1967).  In  prac¬ 
tice,  of  course,  various  limits  are  considered  that  reduce  the  problem  to 
one  of  obtaining  only  a  few  coefficients  of  importance. 

We  employ  here  three  approximations.  The  first  consists  of  the 
replacement  of  the  wave  r  nation  (complete  set  of  Maxwell's  equations)  by 
the  Laplace  equation.  This  can  be  done  if  the  wavelength  is  much  larger 
than  the  dimension  of  the  airplane,  valid  e.g.  for  frequencies  below  500  kHz. 
(This  magnetostatic  approach  was  suggested  by  C.  E.  Baum.)  The  site  error 
is  computed  by  superimposing  the  surface  fields  due  to  orthogonal  incident 
field  components. 

The  second  approximation  consists  of  replacing  the  aircraft  by  a  sim¬ 
ple  geometric  form.  The  most  flexible  3-dimensional  geometric  form  that 
can  be  expressed  by  one  coordinate  of  a  3-parameter  coordinate  system  is 
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the  3-axial  ellipsoid.  Here  solutions  are  expressed  in  elliptical  inte¬ 
grals,  that  need  special  tables  or  computer  programs  for  numerical  evalua¬ 
tion.  However,  if  we  reduce  the  3-axial  ellipsoid  to  a  2-axial  (rotation- 
ally  syiranetric)  prolate  spheroid  or  to  an  elliptical  cylinder,  simpler 
solutions  expressed  in  terms  of  trigonometric  functions  (arc  tangent)  are 
obtained.  In  many  cases  these  simplified  mathematical  models  are  sufficient 
to  study  conditions  on  the  edges  of  the  wings  and  on  the  tail  fin  (using  the 
elliptical  cylinder)  and  the  condition  on  the  fuselage  and  on  the  nose 
(using  the  prolate  spheroid).  Even  if  we  reduce  the  elliptical  cylinder  to 
a  circular  cylinder  the  differences  in  the  surface  fields  due  to  axial  and 
transverse  incident  fields  (the  cause  of  the  site  error)  is  still  signifi¬ 
cant,  as  will  be  shown.  However,  if  we  in  turn  reduce  the  prolate  spheroid 
to  a  sphere  these  essential  differences  vanish  because  of  the  extreme  sym¬ 
metry  of  the  sphere.  Therefore,  the  sphere  cannot  be  used  to  model  the  mis¬ 
direction  effect  of  an  airplane  body. 

Our  third  approximation  is  the  assumption  that  the  aircraft  is  a  good 
conductor,  so  that  the  magnetic  field  is  excluded  from  the  interior  of  the 
aircraft  (by  the  induced  skin  currents).  This  is  justifiable  based  on  the 
fact  that  the  skin  depth  (given  by  6. 6/ AT  cm  in  copper,  where  v  is  the  fre¬ 
quency  (Stratton,  1941))  is  only  0.03  cm  at  a  frequency  of  50  kHz,  and  is 
thus  less  than  typical  aircraft  skin  thicknesses  (mm). 

In  the  magnetostatic  limit,  the  boundary  condition  at  the  aircraft 
surface  is  that  the  magnetic  field  be  tangent  to  the  surface  (zero  normal 
gradient).  (Note  that  these  are  also  the  conditions  in  fluid  "potential 
flow,"  so  that  the  results  of  "potential  flow"  theory  are  applicable.  See 
e.g.  Hess,  1962.) 

We  adopt  a  cartesian  coordinate  system  and  show  in  Fig.  B-l  an  ellipsoid 
projected  onto  the  x-y  plane.  It  is  aligned  with  the  3  axes,  and  has  semi¬ 
axis  lengths  a,  b,  and  c  along  the  x,  y  and  z  axes,  respectively.  (The  z 
axis  is  not  shown.)  We  assume  that  the  top  of  the  airplane  points  in  the  +z 
direction,  and  that  the  incident  magnetic  field  is  parallel  to  the  x-y  plane. 
The  incident  wave  propagates  with  angle  of  attack  eQ  with  respect  to  the 
y-axis  as  shown,  and  with  the  magnetic  field  lines  perpendicular  to  this 
direction.  We  will  consider  the  magnetic  field  distortions  at  several 
points,  such  as  A,  B,  C,  and  D  in  the  figure.  These  can  represent  for 
example  positions  along  the  centerline  on  top  of  a  fuselage  or  vertical 
tail  fin,  or  along  a  wing. 


Av.  - 
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Later  we  will  consider  axially  symmetric  cases  in  which  the  y-axis 
is  the  axis  of  rotational  symmetry.  Thus  we  consider  the  solution  to  be  a 
superposition  of  two  primary  solutions:  One  is  the  "transverse"  case  where 
3  -C  ,  the  incident  field  being  in  the  x-di  recti  on,  perpendicular  to  the 
axis  of  symmetry.  The  other  is  the  "axial"  case  where  9o=S0°,  in  which 
the  incident  field  is  in  the  y-directicn,  parallel  to  the  axis  of  symmetry. 

Before  considering  the  general  solutions  for  the  spheroid  or  ellip¬ 
soid,  let  us  first  consider  the  limit  in  which  b  becomes  infinite.  The 
ellipsoid  then  elongates  into  a  long  cylinder  or  wire  parallel  to  the 
y-axis  and  of  constant  elliptic  cross-section,  as  shown  in  Fig.  B-2. 


A.  CIRCULAR  CYLINDER  (c=a) 

First  we  treat  the  circular  cylinder,  where  c=a  and  b=”,  in  the  axial 
case  where  0o=9O°.  Here  the  field  is  parallel  to  the  cylinder  axis,  and  the 
solution  becomes  trivial.  The  field  is  excluded  from  the  interior  by  a  sole- 
noidal  sheet  current  in  the  surface,  in  the  azimuthal  direction  about  the 
axis.  This  current  produces  an  internal  field  that  cancels  out  the  inci¬ 
dent  field,  but  the  current  produces  no  external  field.  Hence  there  is 

no  distortion  in  the  axial  case,  0  =90°. 

o 

Next  we  treat  the  circular  cylinder  transverse  case,  where  0^=0.  Figure 
B-3a  (App.  B)  shows  the  field  ("flow")  lines  in  this  case.  The  solution  of 
Laplace's  equation  that  satisfies  the  boundary  conditions  at  infinity  and 
on  the  surface  is  the  potential 

4>  =  -  BQ •  ( r  +  a2/r)-cost|>  (1) 


where  r  denotes  the  cylindrical  radial  coordinate,  a  denotes  the  cylinder 
radius,  BQ  denotes  the  magnetic  field  magnitude  at  infinity  (large  r), 
and  ip  is  the  azimuthal  angle  about  the  cylinder  axis  ( 4>=0  along  the  direc¬ 
tion  of  B  ).  It  may  be  readily  verified  that  the  potential  (1)  satisfies 
Laplace's  equation, 


lit  + 1  it  +  J_  lit.  =  n 
2  r  ar  J77 


(2) 


3r 


dtp 


\  • 
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One  boundary  condition  is  that  as  r  becomes  large,  $  approaches 

$  -  -  B  r  cosijj  =  -  Bqx  (3) 

■f 

representing  the  potential  of  the  constant  field  BQ  along  the  x-direction. 
On  the  surface  the  other  boundary  condition  is  that  the  normal  component 
of  the  gradient  vanishes,  that  is, 

~  fr  =  +  Bo(1  — 2")  C0S|^  ^ 

vanishes  for  all  i|>  when  r=a.  Thus  (1)  satisfies  both  boundary  conditions. 

The  radial  and  azimuthal  field  components,  B^  and  B  ,  are  given  by 
the  components  of  the  gradient. 


Br'-f -V<’  -7)-'=“*  (5> 

B*  '  ‘  rl*  *  ■  Bo°  *  7>'s1n*  l«) 

The  distorted  field  of  interest  is  obtained  by  considering  the  x-component 
Bx  at  r=a,  namely, 

(Bx)  =  Br  cos<l>  -  B^  sini- 

=  Bq  +  B0(2sin2i)  -  l)(a2/r2) 

=  Bq  +  (2sin2^  -  1)B0  (7) 

where  the  first  term  on  the  right-hand  side  of  (7)  represents  the  undis¬ 
turbed  field,  and  the  second  term  represents  the  distortion  caused  by  the 
cyl inder. 

Thus,  the  perturbation  field  (addition  to  the  incident  field)  goes 
from  +BQ  on  the  "top"  of  the  cylinder  (where  the  dot  in  Fig.  8-3  represents 
the  probable  position  of  the  detector),  to  -BQ  on  the  "side"  looking  into 
the  magnetic  field.  Assuming  the  crossed  loop  to  be  normally  positioned 
(at  the  dot  position)  on  "top",  it  senses  a  total  field  value  2BQ,  that  is, 
enhanced  by  a  factor  2. 
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Now  we  are  in  a  position  to  calculate  the  misdirection  at  angle  of 
attack  6  .  That  is,  or  top  of  the  cylinder  Bx  is  enhanced  by  a  fac.or 
2,  while  By  remains  unaffected.  Hence  we  have 


D  D  *i 

r  =  tane  =  -  ■=•  tane 

Bx  2Box  2  0 


(8) 


as  the  definition  of  9,  and  therefore 

9  -  e0  =  arc  tan  (1  taneQ)  -  9q 


(9) 


is  the  amount  of  misdirection  or  bearing  error  at  angle  of  attack  eQ. 

From  (9)  we  note  that  the  bearing  error  is  zero  (i.e.,  e=90)  when 
9q=0  and  when  eo=90°.  The  maximum  error  is  -19.5°  (or  +19.5°)  occurring 
when  eQ  is  54.7°  (or  its  supplement  125.3°).  This  follows  from  the  van¬ 
ishing  of  the  derivative  of  (9)  with  respect  to  9q.  The  variation  of 
bearing  error  with  angle  9Q  given  by  (9)  is  shown  in  Table  1, 

At  arbitrary  angle  of  attack  8  ,  the  axial  field  component  gives  rise 
to  skin  currents  flowing  azimuthally  around  the  cylinder,  while  the  trans¬ 
verse  field  component  gives  rise  to  longitudinal  currents,  in  the  ±y  direc¬ 
tions  on  the  two  x-sides  of  the  cylinder.  Thus  the  resultant  skin  current 
flow-lines  are  ellipses  whose  planes  are  parallel  to  the  z-axis  and  the 
direction  of  attack  (i.e.,  skewed  with  respect  to  the  cylinder  axis,  as 
shown  in  Fig.  B-2) . 


B.  ELLIPTIC  CYLINDER  (c*a) 

In  the  more  general  elliptic  cylinder  case,  the  two  principal  fields 
are  again  in  the  x-direction  (eo=0°,  transverse),  and  the  y-direction 
(9q=90o,  axial).  Again  as  in  the  circular  cylinder,  the  axial  field  gives 
rise  to  no  perturbation.  In  the  transverse  case,  however,  the  field  flow 
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is  perturbed  as  shown  in  Fig.  3-3,  where  in  Figs.  B-3b  and  B-3c,  the  flows 
are  topologically  similar  to  the  circular  case.  Fig.  B-3a.  Figure  B-3b 
(c<a)  can  represent  a  flat  wing,  with  the  detector  (large  dot  in  the  middle) 
sensing  relatively  little  field  distortion.  Figure  B-3c  (c>a)  can  represent 
a  thin  tail  fin,  where  the  detector  (large  dot  on  the  edge)  senses  a  large 
distortion. 

It  can  be  shown  that  for  all  cases  in  Fig.  B-3  the  misdirection  is 
given  by 


e  -  eo  =  arc  tan  (R  taneQ)  -  eQ 
where  the  "ratio-factor"  R  is  defined  by 


(10) 


Thus,  R  reduces  to  1/2  when  c=a  as  in  Fig.  B-3a,  and  (10)  reduces  to  (9) 
for  the  circular  cylinder.  At  the  middle  of  a  flat  wing  (Fig.  B-3b)  where 
c<<a,  R  reduces  to  approximately  unity  so  that  the  distortion  is  minimal. 

On  the  edge  of  a  thin  tail  fin  (Fig.  B-3c)  where  c»a,  R  becomes  small  and 
9  becomes  zero  (due  to  the  dominance  of  B  which  becomes  large),  independent 
of  v 


C.  PROLATE  SPHEROID 

A  prolate  spheroid  is  a  reasonable  model  for  a  fuselage  of  finite 
length.  In  this  case  (referring  to  Fig.  B-l ,  where  y  is  the  axis  of  rota¬ 
tion),  the  radii  are  related  by  a=c,  as  in  the  circular  cylinder,  but  b 
is  now  finite.  Hence  the  results  depend  on  the  ratio  t=a/b,  which  is  less 
than  unity  for  a  prolate  spheroid  (but  greater  than  unity  for  an  w  late 
spheroid) . 


—  v .  ^ 
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It  can  be  shown  that  the  ratio-factor  R  in  (10)  and  (11)  is  now 
replaced  by  GR  (see  Eq.  (15)  below),  where  G  is  defined  by 


n  +  \  / 12  2, 1/2  ,,,2  2  ,  .  2  ,2J/2 

G(y,t)  =  (b  -  y  )  / ( b  -  y  +  t  y  )  (12) 

where  t=a/b,  and  y  is  the  y-distance  of  the  surface  point  from  the  mid¬ 
point;  and  R  may  be  defined  by 

R(t)  =  0.5  (1  -  2t2  +  N)/(l  -  N)  (13) 


with  N  defined  by 


In 


1  +  /i-t2 

[rr^: 


Thus  the  misdirection  is  given  by 

9  -  0Q  =  arc  tan  (GR  taneQ)  -  9Q 


(14) 


05) 


where  G  and  R  are  given  by  (12)  and  (13).  Note  that  G  denotes  the  cosine 
of  the  tilt  angle  s  of  the  tangent-plane  with  respect  to  the  horizontal 
plane. 

It  is  of  interest  to  consider  the  limits  t=l  and  t=0.  In  the  former 

2  2  1/2 

(spherical)  limit  R  becomes  1.0,  and  G  becomes  (1-y  /b  )  or  cos3,  where 

3  is  the  angle  of  tilt.  In  the  latter  limit  R  becomes  0.5,  and  G  becomes 
unity,  appropriate  to  the  circular  cylinder  (GR=0.5,  Eqs.  (8)  and  (9)). 

The  variation  of  R(t)  with  t,  and  the  value  of  G  at  selected  values  of  t 
and  y/b,  are  shown  in  Table  2. 

In  Table  3  the  variations  of  bearing  error  with  angle  eQ  are  shown, 
with  t  and  y/b  as  parameters.  Here,  y/b  has  the  values  0,  0.5,  0.9,  0.99, 
and  i.O  (as  in  Table  2),  corresponding  to  the  positions  labelled  A,  B,  C, 
and  D,  respectively,  in  Fig.  B-l  (positions  along  the  top  of  the  fuselage 
where  positions  y/b®0.99  and  1.0  are  both  denoted  by  the  letter  0). 


D.  T-39  MODEL  SITE  ERRORS 


Numerical  solutions  using  a  computer  model  applied  to  realistic 
3-dimensional  aircraft  geometries  are  presently  underway.  The  model 
numerically  solves  the  Laplace  equation,  in  integral  form,  subject  to 
the  boundary  conditions  of  uniform  field  at  infinity  and  zero  normal 
gradient  at  the  aircraft  surface.  The  aircraft  surface  is  approximated 
by  a  multifaceted  surface  where  each  facet  is  a  small  quadrilateral  panel, 
as  illustrated  in  Figs.  B-4  and  B-5. 

In  the  light  of  data  obtained  in  1977  by  the  Air  Force  Flight  Dynamics 
Laboratory  in  flight  tests  of  Stormscope  (see  p.  40)  we  employed  our  com¬ 
puter  to  obtain  a  preliminary  assessment  of  the  possible  influence  of  site 
errors.  The  aircraft  is  a  T-39,  wit;,  the  instrument  installed  near  the 
leading  edge  of  (and  on  the  underside  of)  the  right  wing  tip.  The  model 
portrayed  in  Figs.  B-4  and  B-5  shows  the  wing  modeled  reasonably  realis¬ 
tically,  while  the  fuselage,  whose  detailed  structure  should  be  unimpor¬ 
tant  in  this  case,  is  modeled  crudely. 

The  panels  in  Figs.  B-4  and  B-5  are  labelled  by  letters  A-H,  denoting 
various  sections,  with  A  and  B  on  the  fuselage,  and  C-G  on  the  wing  and  H 
on  the  wingtip.  Each  section  has  12  panels,  with  Nos.  1-6  on  the  upper 
surface  and  Nos.  7-12  denoting  image  positions  on  the  under  surface  (with 
7  under  6,  8  under  5,  ...,  and  12  under  1).  The  Stormscope  instrument 
position  is  on  Panel  G-ll. 


Some  selected  preliminary  results  are  as  follows,  indicating  panel 
location,  maximum  misdirection,  and  angle  of  incidence  at  which  this  occurs. 
For  each  section  we  give  the  optimum  location. 


A  -  9:  +11°  at  135°  (bottom  of  fuselage  at  midwing) 


B: 


D  -  2 
D  -  11: 

E  -  11: 
E  -  2 


(no  good  location,  vertical  plane) 

+16°  at  30°  (top  of  wing,  behind  leading  edge) 

+18°  at  5°  (bottom  of  wing,  behind  leading  edge) 

+8°  at  15°  (top  of  wing,  adjacent  to  leading  edge) 
+12°  at  5°  (bottom  of  wing,  adjacent  to  leading  edge) 

+5°  at  15°  (bottom  of  wing,  adjacent  to  leading  edge) 
+10°  at  160°  (top  of  wing,  adjacent  to  leading  edge) 


'■  w  ■  ■ 
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REAR  VIEW 


(showing  underside  as  well  as  topside) 


FIG.  B-5.  T - 39  COMPUTER  MODEL.  OTHER  VIEWS. 


-88- 


r F  -  5 :  +3°  at  110°  (top  of  wing,  ahead  of  trailing  edge) 

*F  -  10:  +3°  at  25°  (bottom  of  wing,  behind  leading  edge) 

rG  -  10:  +4°  at  80°  (adjacent  to  wingtip,  bottom,  behind  leading  edge) 

*G  -  3:  +6°  at  115°  (adjacent  to  wingtip,  top,  behind  leading  edge) 

G  -  11:  -7°  at  170°  (Stormscope  location,  bottom,  leading  corner  of 

{  0  wingtip) 

G  -  2:  +10°  at  95°  (mirror  of  Stormscope  location,  top  surface) 

H:  (no  good  location,  vertical  plane) 

The  foregoing  represent  optimum  locations  (where  the  misdirections 
are  minimal).  The  misdirections  are  larger  at  other  locations. 

The  following  conclusions  may  be  drawn.  The  optimum  fuselage  location 
is  underneath,  at  midwing  position.  The  optimum  locations  on  the  wing  (and 
in  fact  on  the  whole  airplane)  are  near  the  wingtip  and  away  from  the  fuse¬ 
lage,  either  on  top  and  f.head  of  the  trailing  edge,  or  underneath  and  behind 
the  leading  edge.  The  actual  Stormscope  location  used  was  a  not  unreason¬ 
able  choice  (in  the  absence  of  data  on  site  errors)  but  could  have  been 
improved.  The  reported  bearing  discrepancies  (Baum  and  Seymour,  1980)  are 
consistent  with  the  computed  maximum  misdirection  near  the  wingtip,  of  the 
order  of  10°. 

It  should  be  mentioned  that  these  figures  apply  to  an  incident  field 
lying  entirely  in  the  horizontal  plane.  The  presence  of  a  vertical  field 
component  would  be  associated  with  a  slanted  lightning  channel.  In  this 
way  the  error  due  to  slant  would  be  coupled  with  site  error.  It  is  straight¬ 
forward  to  include  the  vertical  component  in  the  calculations,  i.e.  to  study 
the  site  errors  associated  with  slanted  channels,  in  a  subsequent  extension 
of  the  present  work. 

The  overall  agreement  between  the  analytical  calculations  (e.g.  using 
cylinder  and  ellipsoid  models)  and  the  numbers  obtained  from  the  computer 
model  is  satisfying.  However,  it  is  evident  that  more  detailed  information 
can  be  obtained  with  the  computer  model  calculations.  The  analytical  cal¬ 
culations  are  limited  to  simple  analytical  shapes,  whereas  with  the  computer 
model  one  can  treat  more  realistic  airplane  shapes. 
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Table  1.  MISDIRECTION  ON  CIRCULAR  CYLINDER 


eQ  (degrees) 


Ad  =  misdirection  (Eq.  (9)) 
_ (degrees) _ _ 


0  0. 


5 

-  2.50 

10 

-  4.96 

15 

-  7.37 

20 

-  9.69 

25 

-11.88 

30 

-13.90 

35 

-15.70 

40 

-17.24 

45 

-18.43 

50 

-19.21 

55 

-19.47 

60 

-19.11 

65 

-18.00 

70 

-16.05 

75 

-13.19 

80 

-  9.43 

85 

-  4.92 

90 

0. 

3 


-90- 


Table  2.  MISDIRECTION  CONSTANT  R  VS  t,  AND 
POSITION  FACTOR  G  VS  t  AND  y/b, 

ON  PROLATE  SPHEROID 


Point  A* 

Point  8* 

Point  C* 

t=a/b 

R(t) 

G(y=0) 

G(y=0.5b) 

G(y=0.9b) 

G(y=0.99b) 

0. 

0.5 

1.0 

1.0 

1.0 

1.0 

0.1 

0.5207 

1.0 

0.998 

0.979 

0.819 

0.2 

0.5591 

1.0 

0.993 

0.924 

0.580 

0.3 

0.6054 

0.4 

0.6563 

0.5 

0.7100 

1.0 

0.961 

0.696 

0.274 

0.6 

0.7658 

0.7 

0.8230 

0.8 

0.8812 

0.9 

0.9403 

1.0 

1.0 

1.0+ 

0.866+ 

0.436+ 

0.14V 

★ 

Positions  A,  B, 

C  and  D 

jwi i  in  Fig. 

B-l  (y=0.99b  not  shown). 

+cos  (arc  sin  y/b)  fjr  sphere. 
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Table  3a.  MISDIRECTION  AT  POSITIONS  ON  PROLATE  SPHEROID 
(FIG.  B-l,  t=a/b=0.1) 


(degrees) 

A6(y=0) 

GR=0.521 

Ao(y=0.5b) 

G^-0.520 

A6(y=0.9b) 

GR=0.510 

A6(y=0.99b) 

GR=0.426 

A6(y=b) 

GR=0. 

0 

0. 

0. 

0. 

0. 

0. 

5 

-  2.39 

-  2.40 

-  2.45 

-  2.87 

-  5.0 

10 

-  4.75 

-  4.76 

-  4.86 

-  5.70 

-10.0 

15 

-  7.05 

-  7.07 

-  7.22 

-  8.49 

-15.0 

20 

-  9.26 

-  9.28 

-  9.48 

-11.19 

-20.0 

25 

-11.34 

-11.37 

-11.62 

-13.76 

-25.0 

30 

-13.26 

-13.29 

-13.59 

-16.18 

-30.0 

35 

-14.96 

-14.99 

-15.35 

-18.39 

-35.0 

40 

-16.39 

-16.43 

-16.83 

-20.33 

-40.0 

45 

-17.48 

-17.53 

-17.98 

-21.93 

-45.0 

50 

-13.16 

-18.21 

-18.71 

-23.08 

-50.0 

55 

-18.35 

-18.40 

-18.93 

-23.68 

-55.0 

60 

-17.94 

-17.99 

-18.54 

-23.58 

-60.0 

65 

-16.83 

-15.38 

-17.44 

-22.59 

-65.0 

70 

-14.94 

-14.99 

-15.51 

-20.51 

-70.0 

75 

-12.22 

-12.26 

-12.72 

-17.17 

-75.0 

30 

-  8.70 

-  8.73 

-  9.07 

-12.49 

-80.0 

85 

-  4.53 

-  4.55 

-  4.73 

-  6.61 

-85.0 

90 


0. 


0. 


0. 


0. 


Table  3b 


( 


i 

i 

1 


\ 

I 


i 


eQ  (degrees) 

0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 


40(y*O) 

GR=0.559 

0. 

-  2.20 

-  4.37 

-  6.48 

-  8.50 
-10.39 
-12.11 
-13.62 
-14.87 
-15.79 
-16.33 
-16.40 
-15.93 
-14.83 
-13.07 
-10.61 

-  7.51 

-  3.90 

0. 
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ON  AT  POSITIONS  ON  PROLATE 
[G.  B-l,  t=a/b=0.2) 

SPHEROID 

A0(y=O.5b) 

GR=0.555 

Ae(y=0.9b) 

GR=0.517 

A0(y=O.99b) 
GR=0 . 324 

A0(y=b) 

GR=0. 

0. 

0. 

0. 

0. 

-  2.22 

-  2.41 

-  3.38 

-  5.0 

-  4.41 

-  4.79 

-  6.73 

-10.0 

-  6.54 

-  7.11 

-10.04 

-15.0 

-  8.58 

-  9.34 

-13.27 

-20.0 

-10.49 

-11.45 

-16.41 

-25.0 

-12.23 

-13.38 

-19.40 

-30.0 

-13.76 

-15.10 

-22.22 

-35.0 

-15.03 

-16.55 

-24.79 

-40.0 

-15.97 

-17.66 

-27.05 

-45.0 

-16.52 

-18.36 

-23.89 

-50.0 

-16.60 

-18.56 

-30.17 

-55.0 

-16.13 

-18.16 

-30.70 

-60.0 

-15.04 

-17.05 

-30.21 

-65.0 

-13.26 

-15.15 

-28.33 

-70.0 

-10.77 

-12.40 

-24.59 

-75.0 

-  7.63 

-  8.83 

-18.56 

-80.0 

-  3.96 

-  4.60 

-10.11 

-85.0 

0 


0. 


0 


-'a 
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Table  3c.  MISDIRECTION  AT  POSITIONS  ON  PROLATE  SPHEROID 
(FIG.  B-l,  t=a/b=0.5) 


(degrees) 

Ae(y=0) 
GR=0 .710 

Ae(y=0.5b) 

GR=0.682 

Ae(y=0.9b) 

GR=0.494 

Ae(y=0.99b) 
GR=0. 1945 

Ae(y=b) 

GR=0. 

0 

0. 

0. 

0. 

0. 

0. 

5 

-  1.45 

-  1.59 

-  2.53 

-  4.03 

-  5.0 

10 

-  2.86 

-  3.14 

-  5.02 

-  8.04 

-10.0 

15 

-  4.23 

-  4.64 

-  7.46 

-12.02 

-15.0 

20 

-  5.51 

-  6.06 

-  9.81 

-15.95 

-20.0 

25 

-  6.68 

-  7.36 

-12.03 

-19.82 

-25.0 

30 

-  7.71 

-  8.51 

-14.08 

-23.59 

-30.0 

35 

-  8.57 

-  9.47 

-15.92 

-27.24 

-35.0 

40 

-  9.22 

-10.22 

-17.49 

-30.73 

-40.0 

45 

-  9.63 

-10.71 

-18.71 

-33.99 

-45.0 

50 

-  9.76 

-10.90 

-19.51 

-36.95 

-50.0 

55 

-  9.60 

-10.75 

-19.80 

-39.48 

-55.0 

60 

-  9.12 

-10.25 

-19.45 

-41.38 

-60.0 

65 

-  8.30 

-  9.36 

-18.35 

-42.36 

-65.0 

70 

-  7.14 

-  8.09 

-16.38 

-41.88 

-70.0 

75 

-  5.68 

-  6.45 

-13.48 

-39.02 

-75.0 

80 

-  3.95 

-  4.50 

-  9.64 

-32.19 

-80.0 

85 

-  2.02 

-  2.31 

-  5.04 

-19.22 

-85.0 

90 


0. 


0. 


0 


0. 
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Table  3d.  MISDIRECTION  AT  POSITIONS  ON  PROLATE  SPHEROID 
(FIG.  B-l ,  t=a/b=1.0  for  sphere) 


(degrees) 

A0(y=O) 
GR=1 .0 

A0(y=O.5b) 

GR=0.866 

A@(y=0.9b) 

GR=0.436 

A@(y=0.99b) 
GR=0. 141 

A0(y=b) 

GR=0 

0 

0. 

0. 

0. 

0. 

0. 

5 

0. 

-  0.67 

-  2.82 

-  4.29 

-  5.0 

10 

0. 

-  1.32 

-  5.60 

-  8.58 

-10.0 

15 

0. 

-  1 .94 

-  8.34 

-12.84 

-15.0 

20 

0. 

-  2.50 

-10.99 

-17.06 

-20.0 

25 

0. 

-  3.01 

-13.51 

-21.24 

-25.0 

30 

0. 

-  3.43 

-15.87 

-25.34 

-30.0 

35 

0. 

-  3.77 

-18.03 

-29.36 

-35.0 

40 

0. 

-  3.99 

-19.91 

-33.25 

-40.0 

45 

0. 

-  4.11 

-21.45 

-36.97 

-45.0 

50 

0. 

-  4.10 

-22.55 

-40.46 

-50.0 

55 

0. 

-  3.96 

-23.10 

-43.61 

-55.0 

60 

0. 

-  3.69 

-22.95 

-46.27 

-60.0 

65 

0. 

-  3.30 

-21.93 

-48.17 

-65.0 

70 

0. 

-  2.80 

-19.86 

-48.81 

-70.0 

75 

0. 

-  2.19 

-16.58 

-47.23 

-75.0 

80 

0. 

-  1.51 

-12.02 

-41.34 

-80.0 

85 

0. 

-  0.77 

-  6.35 

-26.81 

-85.0 

90 


0 


0. 


0 


0 
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